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A Idnetic theory approach is used for investigating the structure 
of shoch «ave in a dissociating diatomic gas. An axially elastic 
sphero><»cylinder molecular model vihich permits transfer of rofetional 
and vibrational energy at a collision is proposed for a vlbrationally 
relaxing dantomic gas. The axial elasticity has been Introduced to 
account for the vibrational degrees of freedonu !Phe distribution 
function in such a case depends on a set of 14 variables (g » c ; e fU j 
vdiere B and <j> are reqtiired to describe the vibiational state of a 
mo lee vile f 

The mathematical description used for worlcing out the collision 
dynamics of the axially elastic sphere cylinder molecules, closely 
follows the one used for the rigid spherocylinders by Haight and Lundgre 
The geometry and dynamics of a collision betweai two molecules are studi 
at first, 

!Ihe frequency of binary collisions in the diatomic gas is obtained. 
This is used in writing down the collision integral in the ■Boit 0 mann 
equation. Sbr the BoltzaBon equation, thus fomiilated, an H^Theorem If 
given and through this theorem, the equilibrium solution of the Boltsaaan' 
equation Is obtained, 

Tbr the shock structure- problem, a trimodal and a qUfsirimodal 
ansatz are proposed and a moment method is used to T^nd the wn 



Txbrational energy are used in addition to the conservation eq,nations. 
Ibr Ihe quadramodal approach, the rotationai energy, vibrational energy 

p 

and the C laoment of the Boltamann equation are used to determine the 

weight functions again, in addition to the conservation equations, lUie 

collision integrals in these moment equations are 22 fold out of which 

16 fold integrations are carried out analytically arai the remaining 

integrations were performed numerically using the Monte Carlo technique. 
The set of eqmtions in the unknown weight functions have been solved 
numerically to get the shock profiles for a range of Ifech numbers and 

upstream pressures for nitrogen gas* 

Ibrther, for the case of a dissociating gas the Boltzmann equation 
is modified by including the collision terms representing dissociative 
and recombinative collisions* A separate ^Itzmann equation is writtei 
for the atom species present in the dissociating diatomic gas. dheae 
two Boltzmann equations are solved using a moment method for the shock 
boundary conditions. The trimodal ansatz suggested earlier is used fo' 
the diatomic molecular distribution function and a bimodal distribution 
function for atom specxes. 

In addition to the two species conservation equations, the 
rotational and vibrational energy equations have been used for deter- 
mining the wei^t functions. A^in the collision integrals in these 
moment eqmtions are also 22 fold and have been evaltiated as in the 
earlier case. The shock profiles have been obtained by solving the se' 
of eq^iations in the wei^t functions for various Mach number and upstre 
pressures for nitrogen 
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CHAPTSi 


1 


IKIROIfaCTIO?. 


1 ,1 , This work is concerned with the problem of shock struct-ure 
in a dissociating diatomic gas. In gas dynamics, it is customary 
to assume the shock to be a surface of discontinuity across which 
two uniform states exists These two uniform states are related 
by the Rankine - liigoniot relations. Pfowever, to describe the 
narrow transition zone in the shock y^ave it is necessary to consider 
the effects of -viscosity and thermal conducti-vity. This problem 
was solved within the frame work of the Navi er -Stokes equations 
for a monatomic gas first by Ihylor (1910) and later by Becker 
( 1921 ), Thomas (1944) and others. The thickness of shock vsve 
v;as found to be of the order of a fevi mean free paths of the ^s 
upstream. But the validity of the Navier-Stokes equations for 
the shock phenomenon is questionable as the characteiastic length 
for the shock is so small as to break the continuum hypothesis 
made in the formulation of Ihese equations. Ihe results obtained 
by continuum and kinetic theory approaches for the shock problem, 
as given by Gilbarg and Iholucci (l953)> show roughly equivaloit 
sim)ok thicknesses. I&wever, the merit of a theory should be 
judged by comparison of its results with those given by experiments. 
Observations made for monatomic gases by Talbot (1962 ) and 
Sherman and Talbot (i960), after a careful study of the available 



data, both experimental and theoretical, ehow that tiie ffavrer- 
Stakes equations yield sviffici^tly good results for weak shocks 
(M < 2)i but for strong shocks, (M > z), these results are not 
satisfactory. 

Ob orercome this, BSott-aaith (1951 ) treated the problem for 
the strong shock case, from the molecular point of view by using 
moments of the Bol tzmann equation. He proposed an ansats for the 
distribution function considering the physics of the problem. 'The 
shock region being a thin one, he suggested that molecules from 
both the hot and the cold side might freely interpenetrate into 
this region. He assumed the distribution function to be a weight 
sum of two Maxwellians, one corresponding to the cold side and the 
other to the hot. Sie two wei^t functions were determined 
approximately by satisfying the Boltzmann equation in an average 
sense (moment method). 'Phe results obtained thus were found to b 
satisfactory for a strong shock. Mott-Snith’s analysis was, 
however# restricted to rigid sphere and Sutherland molecules. It 
was further extended to a general inter-molecular inverse-power 
repulsive force law by Muckenfues (i960). 

ibr the case of a diatomic gas with rotational degrees of 
freedom the shock problem was solved by !Ealbot and Scala (l96l) 
from a continuum point of view (modified Navier*Stokes equatinns). 
They defined two temperatures separately one for translational 
degrees and the other for rotational ones and incorpora'ted them 
in the Navier*Stokes equations. For small relaxation times they 
proposed a perturbation technique based on a modified bulk viscosi 



fldiefreae for long relaxation tifflen they used a subsidiaiy relaxatioE 
equation to craaplete the set of equations, “Ilie results obtained b; 
the bulk viscosity approach were found to be unsatisfactory* 
nevertheless j the relaxation model was found to give phenomenologi- 
cally acceptable results which depend on the relaxation time in 
addition to the Mach number* Just as in the monatomic case, in 
this case also, it is reasonable to expect the continuum analysis 
yield results valid for weak shocks only* IJhus Ihe results given 
by the relaxation nod el are limited to weak shock cases only* 

Scala and l^lbot (1963) further extended the two temperatur 
relaxation model to a vibrationally relaxing diatomic gas by 
defining one more temperature to include vibrational relaxation 
effects* Biey also used in this case, two subsidiary relaxation 
equations, instead of one, to complete the set of equations, with 
these modifications the Uavier-Stokes equations were solved along 
with the two relaxation equations for the shock problem. Out of tl 
two relaxation times, they fixed one, the rotational relaxation 
time, and obtained results for different vibrational relaxation 
times and for two Mach nxmabers. In this analysis also physically 
acceptable solutions were obtain si but the limitation of such a 
(continuum) theory, as mentioned earlier, lies in the fact that it 
IS capable of predicting the structure of only weak shock waves 
satisfactorily. 

furcotte and Scholniok (1969) proposed a Beam Continuum Mbdv 
(b*G.M. ) for monatomic strong shocks in hypersonic flows. In this 
the upstream hypersonic flow is taken to be a zero temperature 
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■^ach, the classical desoiip- 

:ontlnum flow is assume 'i 

-’^''11-.. 

Further, this change is supposed to be bro\ight about by 
aollision. per molecule of the beam, The Kavier-Stokes equations 
are used to treat the continuum particles with additional terns to 
account for liie addition of mass, momentum and energy from beam 
particles* Ihe shock profile given by this theoiy is in reasonable 
agreement with some of the experimental profiles given by Oamac 

(1965) and Bobbea and !I!albct (1966)*- - 

Ohis model (b.C.M.) has been later extended by furcotte (1969) 
to the diatomic case (to include the effects of rotational degrees 
of freedom). Bie modified ^Tavier-Stokes equations proposed by 
lalbot and Scala (1961 ) have been used to treat the continuum flow 
in this case. An equivalent haid sphere model was used to specify 
the interaction cross section ’('fetch appears in the ejcpressions for 
the transport coefficients. The results were obtained by varying 
the number of collisions required for rotational relaxation, as a 
parameter for two different Mach numbers. They are in good agreffluesi 
with the experimental results given for strong shocks in nitrogen, 
by Sobben and Talbot (1966). 


1,2. Before dealing with the shock problaa in a diatomic gas from 
a molecular point of -view, one bas to develop a suitable kinetic 
theory which takes into account the internal degrees of freedom. 
There are two available approaches to the kinetic theory of diatomic, 
gases >- the semiMjuantum mechanical ahd the classical. 


JLH tUe seoi— quanttBn mecbanical approach | the classical descilj' 
■Hon fbr 'the translational dogreos of freedom Is retained, tjut the 
internal degrees of freedom are handled quantum mechaaioaliy, The 
difficulties in such a trea-tment lie in the qmntum mechanical 
evaluation of angular dependent in*elastio cross sections. This 
theory was given by Chang and Uhlenbeck (l951 ) and later was 

used by Monchik, Yun and Mason (1963) for evaluating the transport 
properties of gases. To enable a treatment Of a wider class of 
problems in gas dynamics, physically pla\isible kinetic models were 
proposed for a modified Boltsanaon equation, by Morse (1964 ), Brau 
(1967) and Hanson and Morse (1967). These are similar in spirit to 
the B.G.K, model (l954) for the monatomic gas case* 

Further, an analysis close to Grad's thirteen moment approxi- 
mation (1949) for the monatomic gas case has been given by Me Oormio. 
(1968) to get a set of kinetic equations for poly-^tomic gases. He 
used Wang Chang, Ilfalenbeck and de Boer equation (1964) (inelastic 
Boltzmann equation) and four additional moments to take care of the 
internal temperature and the internal heat-flux vector to obtain the 
seventeen moment approximation, These equations are used for 
studying the near equilibrium and relaxation type of flow. The 
restJlts obtained are in very close agreement with those of Wang Chan 
et al (1964) analysis. 

In the classical approach a suitable model is chosen for the 
gas molecule possessing internal degrees of freedom. The tsro 
body collision problan. is solved classically and a modified Boltzsaan 
equation formulated. Such a theory was developed by Curtiss (l956, 
1957) for diatomic gases having rotational degrees of freedom, and 


was fiirther xised for e-vaiuating the transport pro parti os for gases 
with loaded sphere and sphero •cylinder laolecvilar laodels by halher 
and father (1963) and Sandler (l966, 1968). 

1,3. (Phe first attsapt to obtain the shoch structure in a 
rotationally relaxing gas from, the molecular point of -view was made 
by Haight and lundgren (l965t 1967, 1968). They used the kinetic 
theory developed by a classical approach, for the strong shock. A 
trioodai distribution frmction, built up using three Maxwellians, 
describing the Iwo uniform states - the upstream and the down stre£ 
and an intermediate state having no rotational excitation - was 
proposed. A moment method was used to determine the three weight 
fxonctions. This analysis is an extension of Mctt-Sraith*s bimodal 
theory (l95l). Shock profiles for three different cases of 
molecular models (rough sphere, loaded sphere and sphere -cylinder) 
were obtained* As a special case of this theory, it ves possible 
to recover SSbtt-anith mona-tomic results* 

Later the trimodal ansats was used by Erau, Simons and Macon 
( 1969 ) in solving the kinetic model equation proposed by Brau (196'" 
vath the shock boundary conditions. They also tried an alternate 
ansatz based on a model developed by Anderson and Efecomber (l965) 
for monatomic gases. 'fVhile the trimodal theory gave good results 
for both low and hi^ upstream Mach numbers agreeing with the 
experimental data of Eobben and Talbot (1966) for shock profiles ar 
shock thickness, the second ansatz yielded reasonably good results 
only for shock thickness. Uie shock profiles predicted by the sef 
ansatz were not in agreement with the experimental shock profiles : 


high upstream Mach numbers. 

Venkatraman and Morse (l969) used the kinetic model go-van by 
Morse (l964)j to study the structure of a shock wave in diatomic 
gases with rotational degrees of freedom. Ihey assumed a suitably 
modified Mott«3mith bimodal distribution. In this ansatz, in additic, 
to the two weight functions, the parameters of the second Miswellian 
were also taken to be functions of the distance, {Qais waa necessary 
in order to satisfy the energy relaxation equation. Hie three 
conservation equations and two moment equations ( and 1/2 me - E^) 
were used to determine the five unknowns in the problem. Comparison 
of density and temperature profiles with the experimental data of 
Robben and Talbot (l966) is quite good, 

a> 

Recently Machperson (l97l) followed the classical approach 
and a itonte Carlo scheme for obtaining the rotational temperature 
profiles for shock waves in oxygen and nitrogen. He also calculated 
the variation of translational temperature and density inside the 
shock. The intemolecular potential used in this analysis is the 
one given by Parker (l959). ®ie shock thicknesses predicted thus, 
compare well with the eicperimental values obtained by Mairone (l967). 

Venkataraman and Morse (l970) extended their earlier analysis 
( 1969 ) to F^edict the effects of vibrational relaxation on the 
shock 'Structure in dia'feomic gas, A kiaetio model, similar to the 
one given by Morse (1964), with additional terns to include vibration 
effects, was also proposed by tlje®. In addition "to the three 
conservation eq-uationa, they vised c^,( sQC - moments of 

the kinetic equation. Shock structures in oxygen, nitrogen and 


carbon dx -oxide were obtained for high Mach nmbers using this 
nod el » 

1*4» In the present work, an axially elastic sphere --cylinder model 
IS proposed for a diatoiaic gae molecule. Ajclal elasticity has been 
introduced to account for the vibrational degrees of freedom. !Ehe 
distribution function in such a case depends onx, CjOfmjB 
and <f» . Ihe last two of the variables are required to describe the 
vibrational state of a molecule. A mathematical description of the 
molecule is required for working out the dynamics of a collision 
between two molecules. She mathematical representation followed for 
the axially elastic vibrating molecule follows closely, the one used 
by Saight and Iiundgren {l965, 1967, 1968) for the rigid non-central 
particles, 

She geometry and dynamics of a collision between two molecules 
are studied at first. iSie collision integrals in the modified Bolt^ 
equation are formulated. The moment equation corresponding to the 
Maxwell transfer equation Is obtained. As special cases of this, ti 
three conservation equations are written down. iChrough an H-theorem 
the corresponding equilibrium distribution is obtained. 

further for investigating the shock structure in a vibrationall 
and rotationally relaxing diatomic gas, the moment method is used* 
ibllovdng Mott-Smith (l351 ) and iiiight (l965) two multi-mod4l ansatz. 
for the distribution function have be^ proposed » the first one a tri 
using three Maxwellians, the second one, a quadrimodal using fo\ir 
Maxwellians. In the trimodal case a molecule inside the shock Is 
assumed to belong either to the upstream eqxiilibrium state, the 
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iownaijroaja . Maxwellian state or an intermediate equilibrim state 

without vibrational excitations. In the quadrimodal case, in 

addition to these three Maxwellians, ar^ther Maxwellian^ which 

corresponds to an intermediate equilibrium state with only translatior 

degrees of freedom fully excited and rotational and vibrational 

degrees of freedom frozen^ is taken. ilhe three conservation equatior 

effectively give one algebraic equation in the unknown weight 

fimictions, Ob determine three weight functions in the trimodai 

formulation, two momont equations - rotational energy and vibrational 

energy - are used. ibr the quadrimodal approach the fourth equation 
2 

used is 0 moment of the Bolt 2 smann equation, The oolllsion integral,- 
in these moment equations are 22 fold, out of which 16 fold Integra- 
tions are carried out analytically and the remaining Integrations 
were done nunerioally using the i^bnte Carlo technique. Effectively 
we obtain a aet of two non-linear, first order, coupled total 
differential equitionsand an algebraic equation in the three unknown 
weight functions for the txSrtao.dal case whereas for the quadrimodal 
case an .additional differential equation is provided by the 
moment equation, These equations are solved numerically on a 
computer and shock profiles are obtained for various Mach numbers and 
upstream pressiires for nitrogen. 

1*5* As the Mach number increases the degree of dissociation in a 
diatomic gas behind the shock wave becomes increasingly larger* 

Ohus the theory proposed above for the shock structure is valid at or. 
moderately hi^ Mach numbers and it is necessary to include the effec 
of dissociation in the analysis of the shock structure at higher Ifech 


numbers. 
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!I!he analysis outlined, above has to be modified in two ways: 
firstly because of the appearance of a second species, the atom, due 
to the onset of molecular dissociation, an additional Boltzmann 
equation has to be introduced governing the distribution function of 
the atom species j secondly, additional collisional terms appear on 
the right hand side of iiie BoltgEiann equation for the niatomic species 
representing dissociative and recombinative collisions corresponding 
to atom - Boleciile collisions. Ilhus the problem reduces to the 
calculation of the stock, structure in a chemiCially reacting mixture 
of two species, atoms and molecules. In this case also, a moment 
method is used to solve the Boltzmann equations for two species. 53ie 
trimodal ansatz suggested earlier is adopted for the diatomic species 
and the Ifott-Staith bimodal distribution is used for the atom species. 
Two spociee conservation eqmtions and tw> moments of the Boltzmann 
equation for diatomic species (rotational energy and vibrational 
energy) are used to determine the unknown weight functions. She 
shock profiles are obtained for a range of upstream Ifeoh numbers and 
upstream pressxnres in nitrogen. 

In Chapter 2 the collision dynamics of two vibrating npleeules 
is solved. In Chapter 3 is given the fomulatlon of Boltzmann equati 
for a diatomic and also the equilibrium solution for this equatic 
fhe trimodal and quadrimodal analysis for the structure of a shock 
wave in a vlbrationaHy relaxing diatomic has been dealt in 
Chapter 4. fhe dissociaiion and recombination rates are written foi 
intemaol ecular and i!iolecule-~atoa collisions in Chapter 5. In this 
chapter, the formixlation for the shock structure problem con^dering 



the dissociation effects is also given. In Chapter 6 is discussed 
the numerical technique (Monte Gaxlo method) used for evaluating the 
multiple integrals in the mcmeut equations along with the error 
estimates* finally in Chapter 7 are given results and conclusions 
of the present investigations. 



OHAPTSR - 2 


COHESION nmsmoB 


2.1. Molectilax Ifcdel i 

Before developing a kinetic theory for a vibrationally, 
rotationally and translationally relaxing gas, a suitable molecvaar 
model has to be postulated* Jeans (1925), proposed a rigid loaded 
sphere node! for rotationally relaxing diatomic gas molecules, later, 
Hlrschf elder, Curtiss and Bird (l954) introduced a rigid sphere cylinder 
model in the kinetic theory of diatomic gases. Ifowever, one of the 
first few attempts in developing a kinetic theory of non -spherical 
molecules, with rotational degrees of freedom by Curtiss (1956), used 
a general class of molecules which are rigid convex bodies of revolution 
The loaded sphere model is a special case of ihis general model and the 
sphere cylinder, in some sense, also belongs to the general class of 
convex bodies of revolution* 

With the available models the effects of rotational relaxation 
oil transport properties in a diatonic gas, have been estimated by 
Sathor (l963) and ^ndler (1966,1968) and on the shock structure by 
fe.ight (1965) in. such a gas, A classical theory of vibrationally 
relaxing diatomic gas is not available, so far In the literature. 

Por this purpose, a molecular model which permits the transfer of 
vibrational energy in addition to the rotational' and translational 
energies at a collision has IX) be chosen at first. 
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In. tills thesis an axially elastic sphere-^ ylinder model is 
proposed for such a gas. The axial elasticity is introduced for 
the central cylindrical portion to allow the molecules vibrational 
degrees of freedom. This model can be regarded as composed of an 
assembly of two rigid spherical atoms connected by a linear elastic 
spring guided in a rigid hollow cylindrical shell as shown in 
Ti-g* 

The vibrational degrees j as indicated earlier are treated classi 
The justification for this is given in Chapter 5. The molecule is 
assumed to vibrate from an initial equilibrium state. The amplitude 
of the axial vibration is taken to be small compared to the eqviilibriUi 
separation distance between the two stems constituting the molecule. 
The molecule is assumed to be rigid in the transverse direction » so as 
not to excite any bending vihrations in that direction. Jbr a given 
homogeneous species of molecules the frequency of vibrations of a 
molecule is the some for all the molecules. However, the two variabl 
used to describe the vibrational state of a molecule are; the rm-v i mi im 
amplitude, B, of vibration meastired from the equilibrium positions of 
two atoms and the phase angle, ,j) , of the simple harmonic vibrations 

of an individual molecule. The amplitude B is a measure of the vibra 
tional energy of a molecule. 

In addition tc the mass and the two geometrical parameters ( a 
and b) required to specify a rigid sphere -cylinder model, the axially 
elastic spherocylinder model also needs a spring constant, s, for its 
complete description. As the analysis is restricted to small amplitui 
of vibration, a linear spring is an adequate representation. The two 
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T^iriables B and, 0 f describing the vibrational state of a noleenle enter 
not only in the analysis of collision dynanics, but also in the 
definition of a distribution function for the gas. The inportant 
physical constants for sone diatoniic nolecules are shoviG in Table 

(2.1.1 ). 


Surface Representation 


The nolecular surface is described in terns of a unit vector 
e along the longitudinal axis of the molecule and a unit vector n 
normal to the surface, such that the location vector r for a point 
on the surface from the centre of mass of the molecule may be written 
as follows (See Tig. 2.2«1.} : 


r = {b + B cos (vt + } e + a n ) for points on the 

r • n = {b + B cos (vt + 4>)} e»n+aj positive e hemisphere 

w w % m 


r = - {b + B cos (vt + <{t)}e + ap 1 for points on’ the 
r * n = - {b + B cos (vt + <&)} e»n + a \ negative e hemisphere 


(2.2.1) 


r»n = a , e<>n = 0 


for points on Ihe cylindidoal 
surf ac e 


laoics of a Two Body Collision 


In deriving the Boltzmann equation, we shall consider a dilute 
gas and neglect the effects of ternery and higher order collisions, 
as in the case of a monatomic gas. Thus we need solve only a binaiy 
collision problem. The collisional contact time for vibrating 
molecules is finite and is appreciably large compared to the period of 
vibration (Stupochonjco, 1966 ) « However, in this analysis we assume 
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TABLE 2.1, 

1 



■ 

Molecule 


'’2 

°2 

Atomic mass nuraber 

1 

14 

16 

Inter-atomic distance re(A) 

0.7416 

1.094 

1.2074 

Dissociation Energy D(eY) 

4.476 

9.T4 

5.11 

Dissociation Tempo ratuire = D/k(°K) - 

Characteidstio Yib ration temperature 

115000 

59580 

= hv/k(“K) 

5910 

5540 

2250 

Characteristic rotation 

temperature (°K) 

85.4 

2.9 

2.1 

a/h 

CO 

■ 

5.1* 



2J^ 


0 . 025 ^ 



* Galcailations are made for a/lJ =s 5 and « l/52 and ttese are 

£1U ^ 

close to recommended values for nitrogen, molecules dy Haight, 
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the collisioias to be instantaneous. At a collision an impulse I 
is transnltted from molecule 1 to molecule 2 in the direction of the 
coiaoon normal n at the point of contact (lig. 2*3*1 )* line of 

action of the impulse does not almys pass through the centres of 
mass of the molecules. !l?husj in addition to the transfer of linear 
impulse, an angular impulse is also transmitted at a collision, 

Further, the simple harmonic oscillations absorb some energy due to 
the impulse, From the conservation of linear momentum we have, 

In 

c’ = a, - — - for molecule 1 

“1 Hli 

In (2,5* 0 

t ^ 

S'9 ” 'r2 molecule 2 

^ ^ niy 

Ihe angular moiientiaji conservation gives us 

« T-1 

w, = w-j « J- • in X r, 

1 1 *^1 

(2.3.2} 

= (Oo + 0 In X 

The energy absorbed by the vibrational degrees of freedom for molecule 
1 and 2 may be written as 


C&V,H.)j = - (I n) • (ej u Bj sin <j>j) 

(hV,E,^2 = Cl 2^ * (©2 ^2 

Finally we have from the conservation of energy, 

TOjC j j J s(2Bj)^> + ^2 0^ + - 

* 1sC 2B2)^} » {iitii cJ 2 + i . Jj . + isC2B^)^ + (AVE) 

^ % ^2 £2 ^ -2 * -^2 ♦ ^ 2*2 ^^.3 

^ '' -Avv-jo/-. 1, a ft' ' 


1 
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Substituting fron Eqns . (2,3*1 •) to (2.3»3*) in Eqn, (2,3*4* )j 
v/e have for the iapiiLso I i 


I 


(“ g.n)/4 (~ 

V ^ fl. 


TIl^ 


(Ji X r0‘ 


2 J 


(r2 X n)-^ 

+ — } 


11 


2 J 


12 


(2.3.5) 


where 


g - Cc^ + U2 ^ ^2 - ®2 ^2 ‘*'2^ 

( 2 . 3.6 ) 

* -1 ^ ^1 ■ ^1 ^ 

Once I is known in terns of the pre-<!Ollision variables and the 
collision geoaetry, the post collision variables can be detemined 
fron the Eqns. (2.3.1.) to (2.3.3. )» ^hus binaiy collision 
problem is completely solved for vibrating molecules. 

In the formulation of BDlt 2 snann equation for dissociating 
diatomic ^s (see Chapter 5)» we need to consider the collisions 
between an atom and a molecule. Ihe post collision variables for 
such a colli pion can be obtained as a special case of the above by 
putting for particle (atom ) 

5^2 — , B2 0 , (^2 0)2 = 0 (2.3.7) 

2 * 


In Chapter 3j a Beltsamann equation for a diatoiaic gas consisting of 
elastic spherocylinder molecules is formulated and an equilibrium 
solution is obtained for this equation. 
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3.1- Introduction - In this chapter a kinetic theory of vibratiomlly 
relaxing diatomic gas is given. As the end result of this thesis is not 
mxnetic theory per se , a less general view is taken in fomulating the 
Boltzrcann equation. This is especially so in treating the vibrational 
degrees of freedom classically* Strictly speaking, the vibrational 
degrees of freedom exhibit quotum behaviour and the vibrational energjr 
of U identical oscillators (diatomic molecules) is given by Hanck’s 
function! 


Evib 


f hv/kT 
hv/kT . 


(3.1.1) 


vdiere u is the frequency of vibration and T is the translational 
temperature of the gas* In the limit kT >>hv the vibrational 
energy approaches its classical value HkT. Actually the energy 

is already close to its classical value, when kT -v hv * Thus, for 

example, at kT/hv = 0,5j we find 0,724 HkT; at kT/hy * 1 , we 

fxnd 0*928 KkT; at kT/hv =2, v/e find ® 0*979 llkT, 

The value of hv/k for oxygen is 2230'*K and for nitrogen is 
3340°K, (See Table 2.1.1.), 

Sbr a shock wave in a diatomic gas there exists a narrow region 
(translational shock) in v/hich vibrational excitations are practically 
froaen and at the tail of this region the value of kT/hv >1. 
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^jaiere^orei over a significant portion of the shock, following the 
translational shock, the vibrational degrees behave almost classically. 
Hence in the kinetic theory developed here, all the degrees of freedom 
inoluding the vibrational ones are treated classically, 

fiinatioa ~ 

To describo the dynamic state of an axially elastic sphere-^ 
cylinder molecule, a set of 13 coordinates is required. 

The set consists of i 

X - the physical co-ordinates of the centre of mass of the 
molecule 

0 - the translational velocity of the centre of mass of the molecvile 

with respect to the laboratory co-oi-dinate system 
e - the unit vector along the axis of the molecule to fix the 
orientation 

ij} - the angular velocity of the molecule with respect to the 
laboratory co-ordinate system, 

B - the maxiitium amplitude of vibration 

A - the phase angle of the simple harmonic oscillations of the two 
atoms constituting the molecule* 

Thus for a gas made up of such aKJlecules, the distribution function is a 
function varying in the hyper-space q , t where q denotes the space 
defined by ( x, e, u3,B,i|>). Shs probable number of raoleoulee with 
their phase ooordiDates between q and q +dq at any time t , Is given 

by 


f(q;t) dq 
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or 


£(x,c,e,U),B,tfi;t)dx dc de dcj dE do 

where f is the distribution function describing the state of the 
gas and 


dx ~ dx dy dz ; dc = dc^ dcy do^ 


de = sin 0 dE dc ; dw - da )2 '^^3 • 


(5.2.1 ) 


2hus the number density of the molecules at a point x in the physical 
space and at a given instant t is given by the normalisation relation- 
ship? 


nCx,t) = / f (x,c,s,m,B,i{)jt) dx dc dc dw dD d<j) (5*2.2) 

Hie variation of f with time t, is now studied to obtain a 
governing equation for f . In the absence of collisional effects, 
the time derivatives o » of the I 5 phase coordinates of a molecule 
are functions of the coordinates themselves. After an interval of 
time dt from the instant t, a free molecule initially having q as 
its co-ordinate will occupy a point q + q dt in the hyper-space. Now 
consider a set of molecules in a phase volume dq at time t around 

It* 

a point q . At time t + dt the same molecules will be present in 

a phase volume dq* around the point q + d dt, if they suffer no 

J — nr 

collisions during the interval dt, Ifciwever due to collisional effects, 

a fraction of the tolal number of molecules in dq , fall to reach dq* 

* « 

while some molecules which were not present in dq at time t, manage 
to reach dq* , The total number of particles in dq’ , nt time t+dt, 
is related to the number in dq at t by the following equation. 
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fCa+o dt ;t+-dt)ciq ' = f(a;t)tin + {G-L} dq dt 


(3*2.5) 


where the phase volume dq® at o + q dt, t + dt is reifithd . to dq 
at (q , t) by 


dq ' 


* Jt, 


g dt 


dq 


(3.2.4) 


where J is the Jacobian. G dq dt is the number of molecules reaching 

0 *r 

dq’ duiring dt from outside dq j or the number of molecules ’gained’ 

a# 

by do’ due to collie ionsj and h dqdt is the number of molecules in 

dq failing to reach dq» during dt, or the number of molecules 'lost' 

^ -i 

due to collisions. Thus, we have (by expanding the left hand side 
and canoelling the phase volume term on either side) 


= G - L (5.2*5) 

3 qp 

Ihe following are the equations of motion of a free molecule* 


X = c 

C = F - the force per unit mass (3.2,6) 

e = w X e 

^ *T 

J*m=0 

^ •.Mr ^ 

Prom the last equation we have 



J = 



i— ) e e 

^ " 


where 
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so that 


J/j—J-i 


(j) = _ — - — (o) • e) w e 

- „ 

Substituting the above results for in Eqn« (3.2.5) we get 

3f , 9? , r S ^ n 9 

^+c •^+F*^ + — { — ; — (tu*eja! X e £} + — 
9t - 3x - 3c Sd) J, ~ - a; 

^ « J. m 


(3.2.7) 


{(0 X e£} 


+ + HL = G - L. a J . 

3 B Stj) 


(3.2,8) 


Ihe quantities G- and L are some fimctionals of the distribution 
function whose formulation is dealt with, in the following sections, 

3*3. Gollisioa G-eometry and, the Interaction Surface: 


Consider a collision between molecule 1 described by the 

coordinates q and molecule 2 whose coordinates are q_ * Hie 
-1 ~2 

location vectors of the two molecules and ^ differ negligibly 
at the instant of collision. Ibr the type of the geometry of the 
molecule assumed, there exists a common tangent at the point of impact 
of the two collision partners and also a common normal to the surface 
of two molecules at this point, Hcus a collision between two molecules 
1,2 is completely described by the coordinates ( o , e^, ^ # ir]^) 

for molecule 1 , ( f j ) for molecule 2 and the common 

normal h at the point of impact, In ’interaction surface' enclosing 
an ’excluded volume' is defined in the case of a bimry collision, as 
the locus traced by the centre of mass of molecule 2, as it moves in 
contact 7h.th molecule 1, fixing the orientation of both the molecules, 

A point on the interaction sutrface is separated by a vector R from 
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the eeatre of mass of moiecule 1. Bfom the geometry we have 
(Slg. 5.3.10 


5 - 12 > - ^^2 


(30.1) 


go that 

n . dR = n , dr^ - ‘^•^'£2 ~ ^ (3-3*2) 

Hence n is g^so normal to the interaction surface, 5br 1he case 
of a rigid apherocylinder it has already been established by Ifeight 
( 1965 ) that an element dS on the interaction surface depends on 
e.,eQ f n for a given hind of molectsle and is related to an element 
of surface djj on a twit sphere by 

dS * K (Oj, © 2 , n) dn (30*3) 

Extending this for the axially elastic sphero—cylinder case, we 

dS = KCeJ,e2,n,^2»B2,<^J,<f'2) *} (3*3*4) 

where 

2 2 
6 (Cj.n) (e^xn) + (b2 ^2 ( 52 ^!}) ^ 

+ {(bj+B^cos ♦ Cb 2 +B 2 cos <}!2)^ C'S2*B) (5*3.5) 

In Eqn. (5*3-1) and rj functions of ( e^ , n , ) ajad 

(e, , n f B f respectively, where e, and depent^ or " 

rotating molecules. Thus the interaction surface S is variable with 

time. In obtaining the collision frequency we need the normal velocity 
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of a point on the interaction surface. ]?roiii Eqns. (2.2*1 ) and (3«3»1) 
we can write 


' " dt 


(3*3*6) 


where c is the noimal component of the velocity of a point on the 

n,S 

interaction surface S. 

She interaction surface and the normal velocity of a point on 
the surface for a collision between an atom and a molecule may be 
obtained as a special case of the above results by putting 


b2 = 0 , a^^*} » “2“^ ^ ~ ^ *^2 “ 


for the second collision partner viz. atom. 

3.4. GQllislon Erequency ; 

Ito write down the Gain and the Loss terms, G and L in Eqn. 

(3,2.8) we need the number of collisions between type 1 and type 2 
molecules taking place during an interval t,t + dt in the region 
2; , X + dx of the physical space. fhe assimiption of molecular chaos 
is made in finding an expression for collision frequency. Thus the two 

sets of molecules considered in the phase space regions q , q + dq 

-•A ** 1 

and q , q + dq are assumed to be distributed at random without any 
correlations between linear velocity, position, angular velocity, 
orientation, maxirtaim amplitude of vibration and the phase angle of 
vibrations. In other words, all the pre-collision variables of the two 
sets of molecules, are uncorrelated. Ibr long sphero cylinder molecules, 
chattering collisions, in which the pre-collision variables are strongly 
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dependent, are poesidle* Hie effects of such chattering collisions 
jna^ be reduced by choosing a small b/a, 

A molecule of type 2 has a velocity c. ) relative to the 
centre of mass of molecule of type 1 . 5\irther the normal component 
of the velocity of a point, on the interaction surface with respect to 
the centre of mass of molecule 1 is non-aero and is given by Eqn, (3.5»6)* 
Hence the velocity of the centre of mass of molecule 2 in a direction 
normal to the interaction surface and relative to the surface is given 
by 

c-^^c = Cf • n) (3.4,1 ) 

where g is given by Eqn. (2.3.6). 

A collision between molecule 2 and molecule 1 (or collision of 
type 1-2) would take place when the centre of mass of molecule 2 
occupies a point on the interaction surface defined, IChus, the 
molecules of type 2 whose centres of mass are lying inside a cylinder of 
base area dS on the interaction surface S and length - (g n) dt 
noimal to the surface (see lig. 3»3*1 ) would encounter a molecule of 
type 1 during a time interval dt . Hie probable number of type 2 
molecules lying inside this cylindrical volume (- g . n dt ds) in 
physical space and having co-ordinates between q and q + d q is 

^ ( ^2 ’ 5^2 ^ ^2 ■* ^2 * ^2 ^ 6t dS} dcj? ^^2 ^^2 

This also gives the nmbor of type 1-2 collisions tald.ng place diiring 
the time interval dt for one molecule of type 1 present in the phase 
space q , q + d q , 

V X 4S> X 


Thus the total number of 1-2 collisions taking 
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place ia the phase space q , q + d q per uait time is given by 

- 1 -.I - 1 


{ / f dc2 dit)2 dB2 


(3.4.2) 


2!!he number of collisions taking place in a unit volume and unit tljme 
known as collision frequency can be written as 


- / dc2<i«2^®2 ^^2 ^^2 ^1 (3.4#3) 

g.n<0 

Ibr the case of collision between an atom and ,a moleciile the collision 
frequency can be written from Bgn. (3*4*3) as follows* !Eaking the 
particle 2 to be an atcmi we 5aare 0, 62’" 0* iiJ2 ^ ?2 ** 2 • 

Jtirther the area dS^ on the interaction surface S^f in this ease is 
given by putting bg« 0 and ©2 * n 1 in Bgns* (3#3*4) and (3*3*5}. 

Hius the frequency of moleciile-^tom collision is given by 

- / ^^2 d£i de^ dwj dBj d<Jj 

gl.n<0 

where is atomic distribution function and 

li ” 52 ■ ’■il "■ 2i * n ■ Si ®i 

3. 5- Collision Integral j 

A molecule of type 1 loses its identity when it undergoes a 
collision* 2iue the number of molecules of typo 1 lost from the 
phase 'volume q, f q* dq* during a time interval dt due to collisions 
is given by 

{ if £^£^yi~g*n) dS dc2 dw2 de2 dB2 d(J>2) doj dt , (3*5*1 ) 


(3*4.4) 


(3*4*5) 
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vhers yOO - C' for y < 0 (3*5«2) 

= y for y V n 

“Che quantity inside the bracket corresponds to the loss term 
If in the Sqn. ( 3 » 2 * 8 ). ebcosing properly the pre-collision 

variables ( x , c^' , , <’) and ( ^ » '^2 * ^2 *^2 *^2 ^ 

for the two sets of molecules, it is possible to obtain molecules with 
post collision variables as ( x , , a , e. , 

( X , c , fji , Q , * 2 )* such a collision a molecule of type 1 

is generated* Using the expression for the collision frequency 
Eqn. ( 3 * 4 -. 3)5 can write down the total number of 1’- 2’ collisions 
occurring in time dt as 


{/ fJflYC-f' duo de2 db, dr,J 


dt 


(3.5.5} 


itom Iiiouville's theorem of statistical mechanics we have for the 
conservation of phase volume. 


dCjdc^ dwj do52 dS^ 6 . 1 ^ dC'j d^2 " 


'‘ 5 l ^52 ^^2 ^-1 '"“2 ‘-<^2 


(5.5.4) 


It can also be shown ttet 

(g*n) = -Cg’.n) 


( 5 . 5 . 5 ) 


Using these results, the number of 1 2' collisions may be written as 

{/ fi £2 vC-g.n^dS dc 2 dw 2 de 2 8 X 2 doj dt (5.5*6) 

Thi-S also gLves the number of moleculee of type 1 , generated in time 
dt, in space q^, q^ + d q^ * Hence the quantity in the brackets in 
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cozT'espou<is to the gain teiro G in Eqn. ( 5 . 1.7), Ohus 
we have the expression for the gain and logs terms for molecules of 

lying in q^, q^ + dq^ . Krom this we can write the collision 
integral J for the molecules in q , q + dq by dropping the 
subscript 1 and changing subscript 2 to 1 as 


J = / f* 


^ Y("Ghi)d 8 cItj^ (3»5.7) 


where 


^"'$1 f-ri <"--1 ^'^'1 (3,5.8) 

Shis completes the formiilation of the Boltaoaann equation for a gas 
with axially elastic spheix) -cylinder molecules. 

3 . 6 , tfoment Equations ; 

Hie corresponding Maxwell transfer equation for this case may 
be obtained as usual by multiplying both sides of the Boitgmann 
equation (3»2.S) by a molecular property i|j( x , c , e , w , ^ ) 

and integrating with respect to c , w , e , b d • 

•*¥ ^ 

gives us the allowing equation 




f flf. + c In+f? / 


1 


(w.e) (a)>ro)f} 
*» *» 


3 , aAf 9^5-, 


= / J 'vfi dx 


( 3 * 6.1 ) 
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Using the folio vdng notation for the moment of f 


/ dc dfj.' he d3 d(f: = n(x,t) ^(x^t) 


( 3 . 6 . 2 ) 


where g( x , t) is the nmber density of molecules at a point x 
and time t^ the Maxwell transfer equation, for the case where j|, is 
a function of only c r a} and B may be simplified to the form 


3n ^ 
3t 


3 





n {F . 




+ 



) « X e. 



dn d^j de^ dbj d(j>j 


( 5 - 6 . 3 ) 


in v/hioh Bqn. (3.3.4) is used for ds. 2he right hand side of this 
equation, which is the change ( A\{; ^ the property ^ due to 

collisions, may be written in different forms by effecting some changes 
in the variables of integration as done by Haight (1965). fhiis one 
of the useful forms is 


= j / (h' + Tf*! - d) Y('g-n)F dn (3.6.4) 

Ihe primed quantities ( o^' , o^, and (c*, w*, B*, <!>'} 

are the pre-collision collision variables of two colliding moleoiXLes 
giving rise to molecules with poet collision variables ( c^, , (j^p 

and (g,o,B,^). In the present case, a summational invariant ij) 
is definsi, which satisfies the condition 

4- - l|j = 0 (5»6.5) 

Shus, in the Maxwell transfer equation, if is a suimaational 
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invariant:, the collision integral on the right hand side, vanishes. 
A oarefuL exaroination of the collision dynamics, enables us to find 
the follovsdng summational invariants 

1 


= 


n c 
1 2 


c 


1 T ? 

2 



any fionction of e 


any function of m 
o.ny function of s;i 


n 


(5.6.6 ) 


.%• choosing the first three of these, the familiar conservation 
equations for mass, linear momentisa and energy may be written as 
follows: 


|S- . I- . nu = 0 
St 9x 

p C-^ u — . p + p F 


‘3t 

7 


(5.6.7) 


ci.O 


+ u . — ) = - — ^ + 
2 5t " 5x 3x 


5u 


where u , 1, q and p are defined as the Sallowing moments of the 
distribution function f : 

u = i / c £ dr 

-S' 


and 


7 

T 


/ {| m + I Jj mf - ^ sC2l)-} f dT 


(3.6.8) 


= / m c c f dr 


n = / c r c' + i J, mf + ^ 5C2r}^}£dT 
^ “2 - 2JX.2 
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Out of these moments, u and T refer, respectively, to the 
macroscopic velocity and temperature } p and q are related to 
the pressure tensor P and the heat flux vector Q by the following 
relationships : 

_ 2 

q = Q +p (3.6,9) 

and 

t! = P + p U U 

Zt- «* *• 

It should be noted that, though the angular momentum, is conserved 
during a collision, it is not a summational invariant due to the 
two distinct location vectors chosen for the tvro collisional partners. 

3 . 7 . H-fheorem and the Bquilibrium Distribution : 

In this section, an equilibrium solution is obtained for the 
Boltzmann Equation (3*2.8), through an H-Theorem. Zhe rnaaber of 
molecules with dynamic variables q.j in a phase volume + dq.^ 

is given by f(q., , t) dq. . If the state of a gas changes in a time 
interval dt, this number also changes to (f^+ Af-| ) dq.| . The 
fractional change in the number of molecules accompanying a change of 
the state of the gas is ( dq.j or ( A Inf .j ) dq^ . Such a change 

takes place for molecules with different dynamic variables. Thus, it 
is seen that Inf is an important molecular properly associated with 
the state of a gas. Choosing ij* m inf in the Hfexwell transfer 
Equation (3*6,3 we get 

_ — ff 

. Cnc^THr)= lnC~J-)£f, YC-g.n)K dn dr, dr (3*7*1) 
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IHirther, it can be shown that for the case of an axially 
elastic sphere cylinder particle the right hand side of Eq.n. (3.7*1 ) 
is always negative or at most ecLuai to zero. In the case of rigid 
spheroeylinder molecules a similar reeifLt was used by Plight (1965). 
Thus for a uniform case, we can. write the above equation, as 

ff 

/ In K (in dt ( 3 * 7 * 2 ) 

<0 

where 


H - n TnT / f Inf dr (3.7.3) 

Eqh, (3#7*2) is the statement of the H-ThGorem, Sbllowing the 
arguenent given by Chapman and Cowling (l970)^ it can be concluded 
that H decreases rnDnotordcally to a constant value corresponding to 

J U 

the case s 0 * further, it is also knewn that this corresponds 
dt 

to msKimum entixjpy or statistical equilibrium, 33ius for equilibrium, 
we have from Eqn, (3.7*2) 

Inf^' + Inf' = lnf.j + Inf (3.7.4) 

Comparing this vdth Bqn. (3*6*5) it can easily be seen that Inf 
is a summational invariant or a linear combijoation of the summational 
invariants ( 3 , 6 , 6 ). Thus we can 'crite Inf in the form 

5(2B^ j)+c{i; (3*7*5) 

where ! ^*3 » «4 ^5 Ihaatlties which may be 

related to physioai paramecers like number density, temperature 
and velocity of the ^s. Using the noamialisation relationship^ 



( 3 . 6 . 2 ) £^d the definitions ( 3 . 6 . 8 ) of velocity and teu^eratnre 
the Eqn. (3*7.3) can he written in a standard form as 
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B , m . 3/2. 
4-ir ^2nkT 


W 


> pCoij pc Cd) ( 3 . 7 . 6 ) 


where 


/ FCwi {) dwj^^ = 1 


( 3 * 7 . 7 ) 


and &(<{>) is the distribution function for the phase an^e , 
of the vibrational state of the molecules. In the absence of any 
other inforaationy the distribution g( (|s ) for ^ will be assumed to 
be unilorm. Ilius we shall talcs 


G(^) = 


2ir 


( 5 . 7 . 8 ) 


Further in writing Eqn. (3.7*5) the following approximation 
is employedi 



^Toax^ " 



( 3 . 7 . 9 ) 


This is justifiable as the value of for is 113000^K 

(see Ihble 2.1.1). Thus we have the Jtowellian distribution for 
the case of a vlbrationally relaxing diatomic gas. In the next 
chapter two solutions for the plane shocks: structure problem is 
discussed for a vibrationally relaxing dl atomic gas. 



CHAPTER - 4 


SOEMULAHON OP SHOCK SHIUGTCRE HtOBLEM 


4 .1, latrodiJction j 

Hue to the complexity of the collision terms in the BcltEmann 
equation (especially the gain term) it has not been possible (till now) 
to solve the equation even for problems vd.th simple boundary conditions 
iUo overcome this difficulty some models for the Boltzmann equation were 
proposed by making some reasonable assumptions about the collision proc 
by Bhatuagar, Gross and Krook (1954), Gross and Jackson (1959) for the 
monatomic gas case and Brau (l969), Iforse (l964), Hanson and Sforse (196' 
for the diatomic gas case. Tbr the case of nKjna-tomic gases, the model 
equation was solved numerically by using an iteration scheme by 
Iiiepiaann et al (1962), 

Another approach to solve a gasdynamic problem from a molecular 
point of view, consists in choosing a suitable ansata for the xmknown 
distribution function leaving a few free parameters in it. These free 
parameters are determned by lising as many moments of the Boltzmann 
equation as the niaaber of free parameters in the ansatz. The solution 
thus obtained isuoniy an approximate one and depends stron^y on the fo: 
of the ansata used and the choice of the molecifLar properties used in 
the moment equations, 

4.2 Ansata for the Distribution Punction : 

Mott-anith (1951) proposed a bimodal ansata for the distribution 
fmction inside a plane shock wave in a monatomic gas considering the 
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physics of the problem, Ibllowing the two step analysis for shook 
waves in a diatomic gas given by Talbot and Scala (l96l) within the 
frame work of Havier-Stokes theory, Ifeight (1965) extended the Mott- 
Smith binjodal theory to a diatomic case. proposed a trlmodal 

distribution function for the shock pioblem. Both Sfett-Saith and 
Haight used mment methods to solve for the unknown or the free 
parameters in the ansata. ahe aasata proposed by Haight is capable of 
predicting only the effects of rotational excitation on the structure 
of a shock wave. In the present analysis vibrational effects are also 
included in the formulation* lUrther in Chapter 5 the dissociation 
effects are also considered for the shock stixioture problem* 

A saore general ansatz was proposed by Weitasch (l96l) for the 
distribution function in the case of strong departure from themodynami 
eguilibriun as un integral superposition of Maxwellians with suitable 
wei^ting functions. She integral representation may be replaced by a 
quadrature formula as given by jKfarasimha and heshpande (1968) for the 
distribution function as 

N 

I ( 4 - 2.1 ) 

where is a Itawellian bviilt up using ( u^j, n^ ) and is the 
product of the weighting function chosen in the integral representation 
and the weights in the appropriate qmdraturs formula. It can easily 
be shown that the Jubtt-Qaith bimodal ansatz is a partfoiiLar case of the 
general form (4.2.1). Iktendiag a similar representation to the 
diatomic case, we have feights'^ trlmodal ansata as a particular case of 
Hqn. ( 4 . 2 . 1 ). The choice of the parameters used in building up the 
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Ha:xwellians is arbitrary in the general fom of the distribution. 
However, if the choice is made on. the basis of the physics of the 
problem at hand, one can expect a better representation of f. 
further the representation of f may be ijaproved by choosing as large 
an II as possible for a given problem. The value of the number H is 
limited by the magnitude of computational vork involved and the ability 
to guess the parameters, to btiild up the Maxwellians, for a given 
physical situation. 

4 .5 Quadrimodal %po thesis : 

Scala and Talbot (1963) obtained the structure of a shock wave 
in a vibrationally and rotationally relaxing diatomic gas from a 
continuum approach. For this they extended their earlier analysis 
(l96l) by defining a third temperature for the vibiational degrees 
of freedom, Ifowever for the present analysis the experimentally 
established fact ty » is used to visualise a three step 

shock. In the initial narrow region (translational shock) only the 
translational degrees of freedom are assumed to be excited fully and 
in the central aone (rotational shock) the rotational degrees of freedom 
are expected to attain- equilibrium. In both of these regions the 
vibrational excitation is taken to bo negligible, finally, in the 
broad tail of the shock (vibrational shock) the vibrational excitations 
occur completely. TO.th such a picture for a shock wave, it is possible 
to choose four sets of parameters to build up four Maxwellians. (in 
an actual shook these three distinct regions do not exist). These 
Maxwellians are used in (4»2.1 ) to arrive at a quadrimodal distribution 
function for the shock probltaa. 



{Qie foiir Maxwelliaas = '*>4, thus chosen correspond, to the 

following four states of the gas: 

1 , unifom free stream equilibiium condition . 

2* intermediate state existing at the end of the translational shock. 

3, intermediate state existing at the end of the rotational shock, 

4, uniform down stream equllibriiaa conditions. 

(The parameters used in building up these four Maxwellians are obtained 
in the next section, 

4.4* Parameters in the MaxwelllanB : 

Maxwellian f^- ®or a given upstream Ifech nimiber temperature 
and pressure , the parameters for the Maxwellian f.j may be written 
down easily. As the rotational degrees of freedom are excited 
completely, at the upstream temperatures considered (’lljW 400®, 500° K), 

the value of rotational tmperature (f ) for the upstream gas is the sam^:: 
as that of the translational tanperature. However, the vibrational 
excitations of molecules at such low (upstream) touperatures are 
ne^egibie and hence the vibrational temperature T , may be taken to be 
sero, 

OSaus we have 

, I - r. s T « 0 (4.4.1) 

b b ’ b 

As a consequence of sero vibrational temperature, the distribution 
fmction corresponding to vibrational degrees of freedom, in the 
Maxwellian (See Bqn. 3.T,6) becomes a Dirac delta function 
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i.e. lim [{-_!?■ , , 1 ^^^ exD (- = 5 (B) (4.4.2) 

Tv, - 0 »!= Tv. ■ 2k Tv 

1 

further in absence of vibrational excitation, the distribution 
G-((f! ) for the phase angle j> is also dropped. Thus we can v/rite 
f ^ , finally as 


£,=-r-:C 


m 


1 47r^2T;kT. 




2tt}cT. 


-)oXD {- 


r . s2 ^ 2 

inCc - 1 up + JjWj. 
2kTi 


} 6CB) ( 4 . 4 . 5 ) 


Maxwellian f - The translational temperature (T. ), the velocity 
2 tg 

(u^) and the number density (n^) of the gas at the end of the first 
aone may be fixed by using the Rankine-Hugoniot condition for 
monatonic case (y = 5/5). However, the rotational and vibrational 
temperatures (T , T ) of the gas are frozen at v^ues corresponding 
to their free stream values. 


i.e. 


where 



21 

80 

!S. 


(1 


25/7 

'■’i 


•) 







(4.4-4) 


Mf'/Av/S ^Tj) 

Thus we have the Maxwellian f^ as 


(4.4.5) 


— (- 

^ 41: 2i:icT 




.,V2(, 


2irkTi 


■)exp {■ 


2k Tr 




2 

— } 5(B) ?(«.,,) 


(4.4.6) 
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Maxwellian - Pis’ther the translational and rotational teiaperatures 
(l , T ), the number density (n ) and the velocity (vn) of the gas 
at the end of the rotational shock may be fixed using the Eankine- 
Hiiigoniot conditions for a rotation^ly relaxing diatomic gas (y » 7/5), 
while the vibrational temperature continues to remain at its upstream 
value ( aero ) * 

T's 

- = (1 + 

(4.4.7) 


le 


- = ~ Cl - iul - i) 
1 M? ^ 


“1 


u 


Ui 


^ f (1 + 


m: 


Tt = T = T3 ; T « 0, 
''3 3 ^3 


Knowing these (quantities the Maxwellian f^ may be written as 






4Tr 2TrkT, 


27 T]cT ■: 


2k T, 


(4.4.8) 


Maxwellian f, - She ualfoim downstream condition of the gas may be 
determined by using the iiiree familliar conservation equations 


ni = n4 U4 


kT, 


’'iC^ * I h) = "4 f“4 * '^1 

2 -,k„ 

“1 * = “4 * ® 5-^4 


(4.4.9) 


The values of n. ,u, ,T. predicted by these equations are exact for 
perfect gases. However in iiie present work the more realistic values 
for h.jU. and T. for nitrogen, considering real gas effects, oaloulated 
by Athye and Peng (1962) are used* Eeoently haw and ^ristov (1969) 
have given more exact i^lues for the downstream parameters for normal 
shocks in a ta-bular form, for both nitrogen and oxygen. Thus we have 



40 


£ ,=X(_5L_1 3/2 fiL_l (2Bf 


87r2^2TrkT^^ ^2r,k7? ^ 


4 ^^T'a' 


2kT. 


'jS'^^j)(4.4.10^ 


in vifaich all the Ihree temperatures 0). f T and S have the same 




value. 


linally we have the ansata (4,2,1 ) in the form 


f (4.4.11) 

where w^, i- = 1,4 are the unknowa functions of the distance x 
varying through the shock. !Hie boundary conditions on these 

\z> C. 

funotiohs may^ written as 


a= 1 ; WgM « 0 for X -CO 


vn =s w„ = V.' = Oj Vf.*!= 1 for x •♦ +«» 
1 2 5 '4 


(4.4.12) 


Thus the M)lecules in the shock belong to one of the four states 
described by the four Maxwellians f^,l * 1,4. 

4.5 feorosoopic Quantities ; 

The gas dynamic variables like density, velocity and various 
temperatures, Inside the shock may be written in terms of the weight 

functions w. , i = 1 ,4 and the parameters in the four Maxwellians used, 

^ c|ui>T^ 'll' 4 1 ? s 

%■ using the definitions for the macroscopic ^^in Eqns. ( 542 . 2 . ) ahd 

( 5,6.8, ) and the expression for f given by Sqn, (4.4,12) we can 

write down the following expressions; 
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ri. 




j ^ \ 

v/„ + w + w ) 

2 3 4 ^ 


u 


n. 


n 


$ 

Jj 


n T T^ a I 

{ w^ + — w^ + —2. —2. ,- . . .:y t 




7 / u ^ ^2 , 


3 u. 


u. 


+ — )} / (n/a, ) 


a? 

_2 


n^ n T„ n T. 

5= {w^ + T~ W + — ^ ^ W + “i — ^ w } /( ■^) 


n, T, 


(4-5.1 ) 


T n. T, 

^ w^}/(n/ni) 


Thus, once we knov/ the Tariation of the weight functions vdth x we can 
calculate the various (luan titles, using these equations, to obtain 
the shock profiles. 

4.6 fement Eqmtions : 

As indicated earlier, after choosing a form Jbr the distribution 
function for the problem, the unknown parameters (weight functions) 
are determined by satisfying appropriate number of moment of the 
Boltamann equation. The molecular property chosen for the moment 
equation should, of course, be relevant to the problem concerned. 

Ihom the continuity equation (4.1.4), we get a simple algebraic 
equation for w^, i 5= 1,4 as 




( 4 . 6 . 1 ) 
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Bie other two conservation eqxaations (roomentum and energy) along with the 

qnMrintodal ansatz, Sqn, (4»4.12), do not yield any additional independent 

equations for Thus we need three laore suitable moment equations 

to complete the system of eq\iations for Sbr the relaxation flow 

situation, two of the relevant molecular properties chosen are the 

rotational energy and vibrational energy. The third choice closely 

2 

follows the bimodal theory and for simplicity the c moment is lised. 

Sqns. (3.5*3) and (3*5*4)j v/e can ^vrlte down the steady one 
dimensional Maxv/ell transfer equation as 

^ /Qj^i{j£dT=-j/Cijjj+4 -'iJ3^-T|j)YC-g.n) . ffj^Kd^dxjdT (4.6.2) 



(4.6.6) 
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' . , _i 8®’*^ 

SimQarly lor - s ''l r ' 


L 1 2 

A *^1 1 rc t- f— 05 J 

i „uC.T ^ " 

i=l ^ 


(4,6.7) 


and lor ^ 

A 1 ? ( 4 * 6 *®) 

^ Tiu(k T^j) ^ ’ 5 

-onalloea by uslns th® »<=“ "■* 

_ r „ non-^i“ 5 eaQiotial 3 .s j 

53;ig distance • ^ si^here Bjolecule). StirtUer 

/ eauivalent sphere 

ol the upBtrea. g _ of the abore eauatioM 

• vs +vip 1 eft hand side 01 

- number. - -- 

^tteii as functions of si,, t 

uattonscah thus be revfflttss as 

.^ee equations ecu 


^r C ^=7 li 

I Si d^r 2 ijai 
i=i 


(4*6*^) 


1=1 # +v,P free stream 

. r are giveh h® 
where ooeMicteots a 8 

,,, n and the sario^ ratios defined 
MaGh n'ornher -i-i 

4.4) as follows; 


^li 7 Hj 1 

c,i = 'V ' 

2i ^ 1 


(4,6, 10) 




Cti = ^"fi ’ 


and the collisi^*^ 


coefficieu'^ 


are given hy 
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nu u: 


— I . . Co^) 


4 


^2ij 


mn uujXj 




2 . 
to ) 

JL 


(4-6.11 ) 



Eqns, (4.6.9) are 22 fold integrals. The evaluation of these has 
been shown in Appendix A* Out of the 22 fold integrals only 16 fold 
integrations could be performed analytically for most of the cases. 

The remaining 6 fold integrations have been carried out numerically 
using the i^nte Carlo technique^ the details of which are given in 
Chapter 6« Hie set of equations (4.6.1 ) and (4-6.9) have been solved 
numerically and shook profiles have been calculated for four jMach 
numbers (m^ =* 15 1 20, 25, 30 ) and three upstream pressures, (p^ « 

0.1, 0.01, 0-001 atmosphere). The results obtained are discussed in 
Chapter 7. 


4.7 Trimodal hypothesis for the Shock : 

As the relaxation time for the translational and rotational 

degrees of freedom are nearly of the same order (x^ 'V t^) it Is 

possible to further simplify the quadrimodal formulation by a merger 

of the translational and rotatioo&l relaxation zones. This reduces 

the number of Lfexwelllans used in distribution function by one and 

hence w© have a trimodal distribution function. The molecules inside 

the shock belong either to the upstream equilibrium state, f^ , the 

downstream Maxwellian state, f^, or an intermediate equilibrium state 

with no vibrational excitations, f_. The distribution function in thi 

3 


case is given by 
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f = Wjf^ + w^f3 + (4.7.1) 

whexe ^ an<3, are the three unknown weight functions to be 

determined. 

As in the earlier case, a moment method is used to determine the 

weight functions. Ibr the present case, in addition to one algebraic 

eqmtion provided by the three conservation equations given by 

Eqn, (4.6.1) two additional moment equations, viz., vibrational energy 

and rotational energy, are used, ®ias the Ihree equations for the 

present formulation are given by Bqns. (4.6.1) and (4.6.9) (k: = 2,3) 

with Wgss 0 in them. Ihe boundary conditions given by Eqn. (4.4*12) 

hold good for the weight functions w. ,vf, and . Solving these 

13 4 

equations , the shock profiles have been obtained for 6 Mach numbers 
(M^5512,15>18,20,25,30) and three upstream pressure (p^ *0.1 ,0,01 ,0,001 
atmosphere). A discussion of these resialts is given in Chapter 7» 

!Hie next chapter gives a formulation for the shock structure 
problem considering the dissociation effects. 



OHATTSi 


5 


SHX^K SmuCTORE WITH DISSQQIAOIQN EgFBOTS 

5*1* Introdvctlon ; 

We have, so far considered the shock wave structiire In a diatomic 
gas with only the vibrational degrees excited to a limited extent so 
that no molecular dissociation takes place. This postulate is valid a 
only moderately high Ifech numbers. As the Mach number increases the 
degree of dissociation in the diatomic gas behind the shock wave becomei 
increasingly larger. Thus, when Ji| = 20, T^=: 300®K, = 5 mm Hg. , 

the degree of dissociation in the case of nitrogen is approx inately 
'lig» (5.1 ‘I ) gives the degree of dissociation behind a normal shock in 
nitrogen for various upstream Mach numbers and pressures ), It is, 
therefore, necessary to include the effects of dissociation in the 
analysis of the shock structure at these hi^ Mach numbers. 

The analysis outlined in the Chapters 2, 5, 4 has to be modified 
m two wrays. Eirstly, because of the appearance of a second species, 
the atom, due to the onset of molecular dissociation, an additional 
Boltanvann equation has to be introduced governing 1iie distribution fane 
of the atom species. Secondly, additional collision teims appear on 
the ri^t hand side of the Boltzmann equation for the diatomic species 
representing dissociative and recombinative collisions corresponding to 
atom - molecule collisions. 'Thus, ihe problsn reduces to the calcuiata 
of the shock structure in a chemically reacting mixtiure of two species, 
atoms and molecxiLes. ®ie corresponding two Boltzmann equations may be 


written formally as; 


^ T 


if 


fS 


K 
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rt 


= ff.fl + [f.n] - [E.f]dl,-[r..'’],«.,-^tP.F.fjj.oconS.*f'‘>P->=l 


‘recorih 


( 5 . 1 ,’^ 




where f and I* are Ihe distribution functions for the niolecules 
and atoma respectively, ate various collision terms on the right 
hand sides of equations and (5*1.2) are defined as follows: 

i) [f^f] : is the net gain in number of molecules' per uiit tine 

to 'non-reacting' collisions between two molecules 
as already obtained in Eqn, (5*4.5) as J, 

11 ) {f,r‘] : is the not gain in number of molecules per unit vtilme 

per unit tine due to 'non-reacting' collisions 
between a molecule and ^ atom. (This term can be 
obtained as a special case of (i) above. 

iii) Ef.f},.: is the number of molecules lost per unit volume 
' das 

per unit time due to dissociative collisions 
between two molecules. 


iv) number of molecules lost per unit volume 

rlXS 

per unit time due to dissociative collisions 
between a molecule and an atom. 


v) tP»P*^lreccJmb^® number of molecules ^ned per unit volume 

per unit time due to (three body) collisions 
between two atoms In presence of a molecule 
resulting in recombination. 
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vi) is the number of molecules ^Ined per unit yolume 

per unit time due to (three body) collisions 
between two atoms in presence of a third, atom 
resulting in recombination, 

vii) : is the net gain in-* number of *tom per .unit volume per 

unit time due to elastic collisions between two 
atoms, Elis can also be written as a special 
case of (i) above or even directly, 

viii) {f,S’) : is the net gain in amber of '-atoms p*m 

unit volume per unit time due to ’non reacting' 
collisloi^ between a molecule and an atom. !Phle 
can be obtained as a special case of (i), 

ix) (fjf) t is the number of atoms gained per unit volume per 

(tXS 

unit til&B due to dissociative collisions between 
two molecules. 

x) ; is the number of atoms gained per unit volume per 

unit time due to dissociative collisions between 
a molecule and an atom. 

o.^ 

xi) {P.Pif} ; is the number atoms lost per unit volume per unit 

re^coiab 

time due to (three body) collisions between two 
atoms in presence of a molecule 

xii) {J.F.P} ; is the number of atoms lost per unit volume per 

reconh 

unit time due to (three body) collisions between 
two atoms in presence of a third atom, 

She fact that at a dissocia-felve collision one molecule gets con-"'erted into 
two atoms and at a recombinative collision tw5 atoms recombine to form one 
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ffioleciae is used to relate the quantities (iii) to (vl) to qualities 
(ir) to (xii), 

Ihus we have 



1 

J3 

if,f} 

* dis 



1 
' « 

{f|f 

* dis 

(5.1.3) 

[^^^’^Irecomb “ 

h 

{ R,R,f } 

' recomb 


[FjRjR] , = 

' ^ recomo 

i 

{R.R.P } 

^ ’ recomb 



5«2« Dissociation and Recombination Rates ; 


The collisions resulting in dissociation and recombination 
may be represented by the fiillowing chemical dBormulae - 


dis ^ 

+ A + A (5.2.1) 



re comb. 


dis . 

A+A + A (5.2.2) 

^ 

recomb. 


dis ^ 

^ A^A (5.2.5) 



recomb. 

Hie probabilily of sufficiently strongly excited molecule 
diajsociatiiig into two atoms as indicated by forward reaction (5*2.3) 
(photo Ci soci.,tion) is v^ry small ( Zoldovxoh, 196&). iUit; smmv; Is true 
Of the photo- recombination^ represoiitcl by reversed reaction (5.2.3). 
Hence these effects are neglected in this analysis. 




+ A 

2 
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!I3ae rate equations for the chemical reactions mentioned, above may 
be TO'itten as 


dn 


dt 


^2 , dn;^ 

b ‘-1 


n 


A. 


+ k 


r 

1 2 


=4 t— 32 = 




(5.2.4) 


(5-2.5) 


where n^ and are the nvimber densities of atoms and molecules 

respectively and the subscripts 1,2 refer to the reactions in 
Sqn. (5*2*1 ) and Eqn. (5*2.2). Ihe reaction rate constants meiy be 
related by the principle of detailed balancing as - 


•'di ''d, n* - 

« {— r «k(t) 

a S % 

where K(3?) is the equilibrium constant} and n* are n^ equilibrium 
values of the number density of atoms and molecules at the given 
temperature* Qtius, it suffices to Icnow either the dissociation or the 
recombination rate. 


5*2*1. Galculatlon of, the Dissociation Rate ; 

Dissociation is assumed to take place when a molecule collides 
with another moleoifLe or an atom with sufficient energsr, Ohe frequency 
of inter - molecular collisions has already been calculated in Eqn. (3*4 *5 )» 
The atom-molecule collision frequency has also been obtained as a special 
case of the inter-molecule collisions in Eqn, (3*4.4). 

Ihe fraction of collisions in which the total ener^r of the colliding 
partoers exceeds the dissociation energy is given by Zeldovich (l966) as 
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7T ®xp(“ 1^^) (5.2.i) 

o 

y^ere is a nuaiber characterising the degree of participation of 
the internal degrees of freedom in the dissociation process and D is 
the dissociation energy. ^Sie probability of dissociation at an 
energetic collision is known as sterio factor and has to be used to 
write down the dissooiation rate. [Ihus, the dissociation rate for the 
reaction Eqn. (5.2.1) may be written using Eqn. (34.4.5) the above 
two factors as - 


D 


, n So "W 
tf.fJdis ' "o rdw* “ /ffjYt-g.n)K dn dt dTj 


(5*2.7) 


5tor the forward reaction (5.2.2) the dissociation rate may be written 


as 


D 




where (P >3 )? (P^^s^) are the parameters of the molecule-^molecid.e 
^nd molecule « atom collisions, respec lively Experimental values 
of these parameters are available in Stupochenko (l967) for nitrogen 
and oxygen* 


5*2.2. Calculation of the Recombination Rate ; 


Becofflbination is assumed to occur vdien a collision betv^een two 
atoms takes place in the presses of a third particle - either a 

molecule or an atom - which carries part of the bond energy. The 
forward reactions In Bjns. (5.2.1 ) and (5.2.2) correspend to dlesociativ 
collisions while the reversed directions refer to rooombinativo coUls lons. 

I II. T. KAfiPOl? 

3 LSRARY 


i 


,226581 
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recoiabinatioii rates for both the reactions are deterHiiiaed*by 
Sqn, (5*2.6) usihg the piinciple of detailed balancing so that we have 


Cf»P»P3recomb ~ 

(5*2.9) 

fF.'’.«’3Teco* = tf.Flois. 

11103 , mi the collision terns on the right hand sides of the Boltzjnann 
equations (5.‘3.1) and (5.1.2) for the moleciaes and atoms, respectively, 
are determined. 


5*3* Ihe Shocic Structirre Problem ; 

In order to determine the structure of a plane shock wave in a 
dissociating diatomic gas it is now necessary to solve the Boltzmann 
equations (5.1.1) and (5*1.2) subject tc the shock boundary conditions. 
As has been done in the foregoing chapters, a moment metiiod will be used 
to solve the Boltiamann equations, Further, the trimodai ansatz used 
for the non-dissociating ^s earlier, will again be used here for the 
molecular distribution function. Ifott-Smith's bimodal ansatz will be 
used for the atoms. Thus, we vsrite 


f(c) = f j + i?s % * ”4 ^4 

^(£l) “ ^6 * ’’’S ^2 


(5.3.1) 
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in which ) hnd are the upstream and downstreaia 

parsneters for iiie atomic species giTren by 


^ ""l ""P "o, " V " ^ “4 "^4 


C5.3.3) 


^ functions of the distance thiough the shock. 

°2 °2 

5br -file upstream temperatis?es upto 5000*^0 for nitrogen, the degree 
of dissociation is negligibly small so that for all practical purposes 
¥?e Eiay assume =0. Consequently ihe distribution function for the 
atoi^ becomes a Maxwellian 


FCCi) = 2v?5a^n4 ( — exp {- 




(5.3.4) 


with the weight function w^., temperature T and velocity U as the 

b Og ©2 

three unknowns in f. As the masses of a moleciile and an atom are 
comparable the momentum and energy transfer (translational') is efficient 
at a collision., %th this physical picture, the Temperature T and the 

velocity for the atomic species are assumed to be same as that for 

the molecular species, 

i.e, T = T, and U = u 

0^ t o 

vj<. 

where T^ and u are given by tiie Eqns, (4.5*1 ). Thus ^have altogether 
4 unknown weight functions w^, i = 1,3, 4, 5 which have to be determined 
by the moment method. 3he distribution functions must satisfy the tTO 
species conservation equations obtained by using the trimodal ansata 


of Eqn. (5*3.1 ) as 


4 4 


, 1 , "i '‘i 3r “ >= “It ^ i 




2 a4 U 4 ^ = " 2^1 


where 


i6'*3 5> 


W2 0, D / * °ui2 ^ ^i'^23dis 


(5 3.6: 
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'j3ie integrals have all been evaluated In Appendix G, 

'ivvo nope nonents of the Boltsmann equation for the molecules 
are chosen to obtain a complete set of 4 equations for the 4 unknowns, 
ne in the case of non-diesooiating molecules, we shall choose the 
rotational energy and vibrational energy the Ebltamann equation {5.1.1 ). 
Ihe restating moment equations nay be written as ; - 



+ 



HiUi 


Ic 



4 

X C II 

ij=l ^ 



4 

ll I., (y 5(28)") 

i,j=l ^ 3 11 2 


Cy sC2B)2) 


I 'ViWg R (ys(2B)^}3 

i=l i2 ^ 


" i, "rslo, tisC2S)b 

12 


(s 


where I ( ?> ) arc the collision integrals already defined in 
J 

Chapter 4, and 


"o.^CcP) - / ^L4i,?2^'dis 


(5.3.8) 


= 




These integrals are also 22 fold, as in the earlier case and 16 fold 
integration are done analytically in Appendix 0 and the remaining 
integrations are carried out numerically using Monte Carlo technique 
given in Chapter 6, 
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J!jxia the foregoing discussion v^a obtain a set of 4 ordinaryj 
coupledj non-linear^ differential equations Tildcli must be eolyed 
simultaneously for the weight functions subject to the boundary conditions 


•'ffj ^ 1 ; v?5 = W4 = \75 = for x -> -« 

~ ^3 " ® ^ *4 ~ ^OT X 


( 5 . 5 ^ 9 ) 


©le Eqns. ( 5 . 3 * 5 ) and (5. 3*7) are solved numerically for these boundary 
oonditione and shock profiles are obtained for a range of Ufach nimibers and 
three upstream pressures. The results thus obtained are discusBod in 
Gbaptor 7» 



OHmSi - 6 

ro^CAI. EYALUATIQH U^TEGBAT.S 


6,1 Introduction : 

3iie collision integrals in the jooment equations encountered in 
the preceding tno chapters are 22 iold and in most of the cases 16 
fold Integrations were carried out analytically. It ms not possible 
to do the remiziing 6 fold integration analytically and hence a numerical 
approach was resorted to. In any of the standard techniques of multiple 
quadrature, the number of points required to evaluate an integral @:)e 5 
up in geometrical progression as the order of the manifold integral. 

Siis exponential increase never occurs in case of Monte Oarlo evaluation. 
Ihus for the numerical integration of manifold integrals this method is 
very economical. Ihe Mcnte Carlo technique gives a reasonable estimate 
of the integral and errors in the integration can also be estimated during 
the process of integration, 

6*2, Monte Oarlo Integration : 

A typical integral which has to be evaluated may be written 
foimally from Appendix A or G as 

I = / dB, (6.2.1) 

0 0 0 0 0 0 

fbr a ISSbnte Carlo evaluation, of this integral the limits of integration 
are changed to 0 to 1 for all the variables by using linear transformations 
for then* Qius we have the integral 
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1 


dx/^dy/^dz/^du/^dv/^ H(x,y ,z ,u,v,w) dw 
0 0 0 0 0 0 


( 6 . 2 , 2 ) 


where 

H(x,y,z,u,v,w) = 4 ^"^ FCB,Ej,i^,(J)j,B,e2) 


(6.2.5) 


The method consists of evaluating the integrand H at N randomly 


chosen points i = 1 ,1T lying vdth-in the domain 

of region of integration and finding average of these F values. This 
v3-lii© itiSGlf is Sin eS'tisi 3 ''t 6 th© xii'bcgf^l* 

i,e. I == leN - ^ H(Xpyi,Z 5 ^»Ui,Vi,w.) (6,2,4) 

Jbr this calculation the points are chosen as random 

numbers belonging to 6 independent series of unifomly distributed random 


numbers between 0 to 1 . Bor the numerical computational purposes the 
random digits available in the Rand Ctorporation Table have been used in 
this thesis. 


6.5 Error Estimation : 

The variance of H in defined by 

dT = h2 - V (6.3.1 ) 

0 

where 

dx = dx dy dz du dv dw, H = j ^ 

0 

and the standard deviation is given by 

= a/t¥ 


(6.3.3) 
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the estimation of Integral Ij then a _+ 15^ error in this error 
estimation allows the overall error in I to tafce a value between 
3,5% and 11.5^* 


In pratical calculations the choice of W depends on the error 
that can be allowed, Lordi and Mates (l970) evaluated 6 fold 
integrals encountered in their analysis, using Monte Carlo technique 
and obtained satisfactory estimates for the multi-dimensional integrals. 
They us©i 400-500 trials for different cases depending on the error 
criterion, Ihree typical plots of the estimates and verstis 

17 are given in Pig. (6.3.1) for the case of three collision integrals 
far , Rotational Ihergy and Vibrational Etaergy. It is found that 
600 trials are required to limit the overall error to _+ 10^ in the 
last two cases* Hcwever, for the moment integral as many as 1000 
trials were necessary to restrict the error to _+ 1C^» In a few cases 
of the o^ moment integrals, the errors were as high as 26% for 1000 
trials • 

Ibr evaluating a typical six fold integral 1^^ (f) or 
given in Appendix A or G it takes about 3 minutes on the IBI 7044 > 
y/ith 600 Monte Carlo trials* Ibr the cases when i "the integrals 

are 7 fold and the same procedure is adopted for evaluating this 
integral. !Ihe computer time taken for these cases is about 3^ minutes. 

As the integrations are carried out approximately the shock profiles 
calculated may also be expected to be approximate to the same extent. 

In the next, chapter a discussion on the results and conclusions are 


given. 



CHAiTER .7 

DISCUSSIONS AND OOHCDUSIONS 


7»t* Shock Structure in a Vibra'tionally Helaxing Oas 

An attempt was made in Chapter 4 to investigate the effects of 
vibrational relaxation on the structure of a shock wave and two 
formulations were presented for diatomic gases. Both the trimodal 
and the quadrimodal theories are extensions of the Mott Smith 'bimodal 
approach and the trimodal theory of Haight for a rotationally relaxing, 
dxatomic gas. Ihe various shock profiles for nitaxigen predicted by 
the trimodai theory of Chapter 4 are shown in Figs. 7*1*1 to 16 and 
the results given by the quadrimodal approach are presented in Jigs. 
7,2.1 to 9. 

7*1*1* Shock Profiles (^Crimodal [Theory): She trimodai calculations 
are carried out for the upstream Mach numbers ranging from 12 to 30 ai 
for the upstream pressures of 0.1, 0.01 and 0.001 atmosphere. Ihe 
density profiles for these Iilach numbers ore given in Tigs. 7*1*1 to 6 
It is seen from these profiles that for a given upstream Each number 
tile shock transition takes place faster with the decrease of the 
upstream pressure, whereas for a given upstream pressure the shock 
becomes thinner with the increase of the upstream Mach number. 

!Ehe corresponding temperature profiles are given in Ilgs* 7.1* 
to 12, fhe dependence of the width of the gone of temperature varia' 
on the free— stream conditions, is same as in the case of density, Dn 
all the cases, the translational temperature initially rises at a fas- 
rate to a maximum and then decreases gradually to its down-stream 
equilibrium value, thus exhibiting a pronounced temperaturo peak. The 
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the rotational profile lags behind the translational one* it shows a peak 
for Maoh numbers abo-ve 20. Hiis overshoot is Tery much less pronounced 
and it occurs slightly later when compared to that of translational 
temperature. the vibrational temperature monotonically increases 

CO its down -stream equilibrtxim. value and further it is found to lag behind 
the other two temperatures. Ibis lag in the vibrational tempera twe is 
significant for Mach numbers above 20* 

Bee above pattern of the three temperature distributions may be 
explained as follows; She translational degrees being the fastest of 
the three, the corresponding temperature attains its equilibrium value 
(monatomio) first with-in a short distance ( 'v Xj^ ). This is also true 
of the rotational degrees of freedom which is slower than the translational 
but faster than the vibrational degrees of freedom. Tinally the slowest 
of the three namely the vibrational degree absorb energy from the faster 
ones and relax to the equilibrium state gradually. The degree of lag of 
vibrational temperature at lower Mach numbers is small mainly due to the 
classical nature of the vibrational mode assumed. 

In ligs. 7.1.13 and 14 are shown the density and temperature 
profiles with the upstream temperature as a parameter. It is seen from 
these profiles that for a given , an increase in the upstream temperature 
reduces the extent of the shock due to higher absolute temperatures inside 

the shock. 

A typical variation of the weight functions w^,w^ and w^ is given 
in Tig. 7.1.15. The functions and show a mono tonic behaviour 
within the sho ck, whereas as expected lias a maximum inside the 
shook with zero values at the two ends, Ibr > 20, w^ is found to 

Ihis fact is reflected in the non-existence of a peak 


take small values. 
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in the rotational temperatiire and a faster approach of the vibrational 
temperature to the rotational temperature. Ehie is perhaps attributable 
to the lower population of the molecules charaoterised by f^ inside the 
shock, .as the dovm<*Btream condition (f^) is attained faster due to the 
classical behaviour of the vibrational degrees of freedom assumed# Ihls 
xs not so at higher llach numbers as this assumption regarding the vibratioi 
degrees of freedom is valid to a better extent. 

A comparision of 'the, theoretical, profiles with the experimental el 
profiles given by Sobben and lalbot (l966) is made in Hg. 7*1»16 for 

12«9* !I3ae density profile predicted by the trimodal theory compare 
very well with the experimental profiles over more than half the distance 

throug^i the shock and differ sli^tly over rest of the distance. But 

.t 

for a shift in the calculated rotational temperature profile, also: compares 
very well vdth the experimental profile. 

7.1.2 Shock Profiles (Qoadrimodal Theory); The density profiles for th< 
quadrimodal case are given in Tigs# 7.2.1 to 5 and the temperature profilse 
in Tigs, 7.2.4 to 6. The calculations were done for 15, 20 and 25 and 
for the upstream pressures indicated earlier. (Bor M^«=30, the shock profile 
ai*e obtained at p^— 0,1 atmosphere and are shown in Tig. 7.2.7.) 
profiles in this case are qualitatively same as in the trimodal case. The 
translational temperature overshoots are some what hi^er in this case 
than for the earlier theory. 

In Tig, 7.2,9 a comparison of the various profiles given by both 
the theories is made (while comparing p w 0,5 is matched for both the 
cases v/here p — (p / (p4 *• Pj )♦ density profiles match closely 

upto a certain distance from the upstream whereas, tl^ temperature profile 
compare fairly closely but for the over shoot in Tig. 7*2,8 shows a 

typical tion of the wei^t functions w^. The fhnetrone and 
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have monotorsic hehaviour as in the earlier case and and possess 

a mepcimuEi vdthin the shockj haid-ng aero values at the two ends. 

7,1.3. Shock Thickaess: The following three shock thiclmesses are 

calculated for the vihrationally relaxing diatomic gas using the profiles 

given by both the theories: density thiclmess ( Ap )» vlhrational tempera- 

,ure thicloaess ( Ay ) and the rotational temperature thickness ( Aj^ )» Ih 

all the cases the maximum slopes of the corresponding profiles are used. 

These thickness indicate the extent of the shock when the actual profiles 

are replaced by linear profiles with the corresponding nmimum slopes. 

The variation of the reciprocal shock thicknesses with the upstream Mach 

number is shown in Tigs. 7.4.1 to 5* Tbr nondimensionalising the shock 

thickness a length L* adopted by Robben and Talbot (l966) Is used instead 

of the free—streairi • 1'^' is seen from the :S.gures that the density 

thickness is lower than the two temperature thicknesses at higher Mach 

the 

numbers , ’This -i a expected due to the presence rotational and vlbrationaj 

relaxiation gones' at higher Mach numbers. The quadrimodal approach predict 

lov/er values for the shook thicknesses. The experimental values of Ap 

and A at M,= 12,9 given by Robben and Talbot (1966) are also shown 
R ■ 

for comparision. The theoretical Aj^ compares well with the experimental 
value at M^=12.9. finally, for a given the shock thicknesses are found 
to decrease with the upstream pressure. 

7.2. Shock Structure in a Dissociating llatomic Gas t 

Hie analysis given in Chapter 4 v/as further extended to include 
the effects of dissociation in Chapter 5. Oie shook profiles thus 
obtained are presented In Pigs. 7-3,1 to 9. 

7.2.1. Shock Profiles with Dissociation Effects! In Pigs. 7.3.1 to 3 
the de«ret of dissooiatioa, total density and ffloleoular density profUes 
are given for !l, 15t 20 ani 25 roapectively. Ihe profiles are also given 
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for M^= 25 wi-th aa parameter in Ifg, 7.3.4. Hhe eorresponding 
temperature profiles are given in Ilgs. 7.3.5 to 8. In all these 
profiles the total density is found to increase monotonioally whereas 
the nolecular density has a maxi mum inside the shock. !Ehe difference 
of these densities is related to the degree of dissociation by a = i , pj^/p 
Ihe study made by Ibrranee (l96l) indicates the possibility of maxima 
for both p and a * However in the present case inspite of a peak in 
the profiles for ^ the profiles for degree of dissociation have a 
mono tonic behaviour, 

Ibr a low a ('v 0.2) the three temperatures are 

found to equilibrate in a narrow region y/here the dissociation is 

practically absent, ISiis region is followed by a very broad clearly 

distinguishable tail in which the degree of dissociation attains its 

the 

equilibrium value. In this broad regioE^changes in p , ’’t' and 
are very gradual whereas the molc-oular density and the degree of dissocia- 
tion change significantly (see Jig, T.3.1, 2, 5 and 6), 

Ibr a high «('b ,9) the two regions lose their identity and 

merge together. However there exists a broad zone in which the throe 
temperatures attain nearly their equilibrium value. Eiis is followed by 
a tail whose width is also of the ssime order. In this region the only 
quantities which vary significantly are the two densities and the degree 
of dissociation, (see Jigs, 7,3.3 and 7). 

finally for a moderately high a ( "vi 0,6 ) there exists a broad 
region in which the three temperatures and the total density attain 
almost their equilibrium value. This is followed by a very broad region 
in which the molecular density and the degree of dissociation attain 
their equilibrium values as in the case of lovir p , In otherwords the 



65 


merger of the two zones vail have just begun, (see Elge. 7.3,4 and 8), 
ihese figures also indicate the effects of the upstream pressure 
on the various profiles. As found in the case of shock structure without 
d.x3sociation, in this case also^ with the decrease of p^ the extent of the 
shock decreases, 

7*2,1. Shock Bmckness with Idlssociation Effectsj She definition of 
shock thickness should have some bearing to the kind of relaxation 
phenomena predominent, Jbr the shock in dissociating gas the degree of 
dissociation of the gas has to be incorporated in the definition of a 
thickness, Jbr the present case in the absence of any other definition 
in the literature, a shook thickness A is defined as a distance from 
the begining of the shock transition to a point inside the shook whore 
the degree of dissociation is 90 ^ of the down-stream equilibrium ralue, 
(a similar definition has been adopted in the c»se of ionising shock 
waves in mono tonic gases). Ihe position where shock transition starts 
is of course not very well defined and to this extent there is certain 
amount of arbitrariness in the value of A chosen. 

Jhr a given upstream pressure the equilibrium degree of dissoci- 
ation increases with Ihe Mach number. JUrther with the increase of 
equilibrium degree of dissociation the shock thickness decreases, liis 
fact is reflected in lig. 7»4,4 which shows the variation of a and 
1 /a with for nitrogen 


7 



7 *5 • Oonc lug ion s ; 

For the comparison purposes hardly any expericiGntal shoclc 
profile data for nitrogen in the Maoh number range considered is 
available. Thus no cSefinice conclusion, can be made as to the validity 
of the theory, Hovfsver the two fonmlations-quadriiuodal and trii-iodal- 
given in Chapter 4 are capable of predicting the shook stnicturo in a 
vibrationally relaxing diatomic gas. Out of these tvro, the quadrimodal 
theory is better than the trimodal one especially for the Hach niunber 
range 15 to 20, Further the analysis presented in Chapter 5 incl-uues the 
dissociation effects in the gas and hence is valid for higher Hach number 
upto 5c*. For > 50, the ionisation and photo-reoombination of the 
atoms and photo-dissociation of the molecules become appreciable. These 
effects are not inclT;uied in iiie present foraialation. 

As particular cases of the quadrimodal theory, for vibrationally 
relaxing diatomic gas given in Chapter 4f it is poesihle to recover, hoth 
the trimodal theory proposed by Haight for rotationally relaxing diatomic 
gases (w^=0) and the Mott Snith himodal theory (w^tsw^aO) for inert gases* 
Further to enable to carryout meaningful calculations for shook 
structure, the parameters b, a and s for axially elastic sphorocirlinder 
models corresponding to real diatomic gas molecules should be available* 
In the present case the values indicated in the Table 2,1.1 wore used. 

The shock structure analysis given in Chapter 5 for the 
dissociating gas may further be extended for a mixture of atoms caid 
molecules of the same -gas. In such a case the upstream atomic number 
density is nonaero and hence the full bdmodal ansata for the distribution 
function for atoms has to be retained* Correspondingly, aPP-~opriat© 
number of relevant moments of the tvTO Soltsmann equations have to bo 
consa derod in Eddition to the equations of Chaptor 5. 
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APPMDIX - A 


COIilSIOM INTEGRALS ^(2 ) 

- ■ - I - -- -- 4.J 


The collision integrals in the moment equations are given by 


lijft) = 


- Y(**g»n)K dn dx dr, (A,1.l) 


where 


g.n - {(Cj^ + X - Ej V sin e^) 

(c + wxr-Bv sin (;> e) }. n 

-V — «* %. 

K dn elemental area on the interaction surface and y{x.) = 0 
for X ■«< 0, “ 1 X > 0. 

Jtam the dynamics of binary collision (Chapter 2)we have for the impulse 
I transmitted at a collision from molecule 1 to molecule 2 


I = - 


g.n 


(A.1,3) 


where 


5 = 1+ -^(b+B cos (}>)^Cexr!)^+Cb+Bj cos Ce^icn) (A, 1.4) 

2i.j 

further the post collision variables are given by the following 
equations 


c'=c+in;c»=c,in 
*1 "1 m - 


m 


Sj = »J * j- ct> . Bj cos tgi « n) 


(ii = a - — (b + D cos 4) (e n) 
- - 


(a.1.5) 
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Tae collision integral given by Eqn. {A,1 ,i) Is a 22 fold 

2 2 

integral. Without making a particular choice for f (c , 1/2 J. m or 

X i 4, 

2 

1/2 s( 2B) ) it is possible to carry out 11 integratiom analytically. 

The dependence of the Integrand on e, e, , n , u and w, is limited to 

•'S 

the following combination of them; 

e»n , e^»n * excj^'n, w.e (oj^j 1 , m^oe^t w.exn 

and Uj • 

Thus the vectors o) and m. may be expressed inteims of e , and n^^ 
as follows (see Tig. facilitating the integrations in w and 

space. In such a case we have 


exn 


exn 


w = w{e cos 6 + ( ; ■ -■ ■■ cos ^ + e X 


! exn I 


exn 


sin iji) sin 0 } 


e^xn 




( A . 1 * b ^ 


21 , Oil {ej cos 01 * Ctt^^ cos ♦ gj x sin ^psln ep 

1 *r i > ^ ± • 


= 


Vi 
dm 

and 


= w 


u sin 9 , sin 

|2 sin 8 do) d6 dij) = dm^ dto^^ d^ 


dMj = Wj sin 9j d6j d\J»j = dujjj^ doj^ d^*^ 


Further representing e and ^ in spherical co-ordinates 


we get 


,,neose* {^cosc *5 gl si» e) sin P 
* - nxgi pxgl 


ej^ « n cos 02 + (s cos + t sin Ej) sin 0 ^ 
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and 


de - sin 0 dB de 


de^ = sin de^^ 


CAa.7) 


where s , t; and n form an orthogonal set. Consequently we have the 
folloY/ing relationships useful for the integration 

e»n = cos B ; 21*2 ~ "■ ^1> 


and. 


exej«n • sin B sin Bj cos e 


w.exn = m sin 6 cos if; sinB 


Wj,»ejXn = sin 9^ cos i|;j sin Pj^ 


(A.1.8) 


Using these relationships we can rewrite Eqns. (a, 1.2) and (A.I.4) 
as 

g.n = v»n + cos 4i) Wj4_ cos ^ sin Bj - (b B cos (J/) cos if* 


sin B - \J sin cos ^ B sin ^ cos B) 


(a. 1.9) 


^ 1 + {(l+Wj COS 4^)^ 3in^Pj^+(l+B cos <p)^ sin^ B) (A.I.IO) 


2J 




V * " G 

^ •f 1 


(A.1.11 ) 


Now the integrations with rpspoot td ? and e-j^ may be 
Carried out^ as the dependence on them is only throu^ K» !Phus 
we have 


/ K de de^ 5 ^ (2ff)^ K 


(A.1.12) 
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where 

B = B/b , = B^b 

and 




K = 4a [1+2 ^Cl+B cos <,)5(cos a3sin^ 6 + (1+B^ cos (f^)5(-cos 


2 2 

sin (1+E cos <J>) (1+F^ cos <J>j^)sin^0 sin^g^dCcos 3) <S(-cos 3)] 

CA.l* 

In Hie integrand, n appears tlirou^ the tena v.n « Eepresenting n 

~ *1 — w 

inteuns r and i we get 

V vxi V v^i 

n = ~ cos a + cos 5 + - X ■,- * ■ sin S} sin a 

V vxi V vxi 


vTi = V COS a 

•# 


(A.t.U) 


and 


dn sin a da d3 

*4 


Ihe dependence of the integrand on 6 if any is only through {|} and 
we can write^ at this state 

2Tr 

/ {i^} dd - 2it {f} (A*1.15} 

0 

Ihus we have the collision integral given by Eqn. (a,1.i) as* 

I .(I) = -(27r)^ / TfT K sin a da dr d7, (A,1,16) 

■i'J **44 ^ 


13 ) 


^ere 
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dx 5= bd c B d3 d^ d (1} 


d-z^ = b dCj^ -^'ii ^1 *^^1 


(A. 1.17) 


and 


-F = !1E_ C— 3 — (jiL-.) %v - 

P r, a 2TikT^ 2TrkT- SrkTv ^ 2^7* 

% ^F '^P 


nr 


J, (1) 


2kT- 21cT^, 

in r^iich F(w^i ) satisfies the fbll owing condition 

CO 

J FCw^j) dwjj » U 


(A.1.18) 


(A.1.19) 


]By using this ( A.1.19) two trivial integrations with respect to 
and (jjj may be carried out straight away, 

She integrations in the velocity space c , C-| W carried 
out by osi ng the stand.ard transformations 


c = u. i + 7 - G 


^4.- 

ti 


IT. - ^ S .+ T, . 

^ " ti t3 


£,= u. i + V S 


(A.1,20) 


such that 

dCdC^ =dVdG- 
^y inspection it can easily be seen 


(A.1.21 ) 

that integrand is dependent 


on Y 


_ only through the exponential te«, ^ua performing the 
inte'gratioae In V epeee .e get the eolUeion integral given hy 

(a,i, 17 ) as 
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!..($) = -C2it3^ L_{ 


m 


(8^2)2 2irk(Tt,+Tt ) 




J- 


-> (* 




27rkT^ 2trkT^ 
H 3 


4 


^J^^l/2 b2 /ITT^YC-g.n) exP {- 


mG 




2k (Tt.+Tt.) 
1 3 


3 ^ s(2B)^ sC2B0‘ 


2kT_ akT^, 2kT« 2kT-. 

»i ^3 


-} K sin a da dr’ dr. 


{A,1.22) 


where 


dr ’ = dG Wjjdw^ di/j 

d7^ » sin P sin pj d!' 3<fi dp dPj 


(a. 1.25) 


and the Jbllowing result is used 


T T.. T 

, f T1 H ^3 ,,2't f2T5k 1 5^3/2 


• t * ' t. 

H ^3 


(A.1 ,24) 


Further transforming ) and Ai ) to a cartesxan 

co-ordinates, after suitably non -dimen sionali sing them we have 


2k T, 


m. dw, di^ - 


‘'“x ‘‘“y 


T’rj 


(A.1. 25) 


WlL dmj^_, = — j- 


dm^^ dwj^y 


lidi.ere 


“x =/iSy“x“’*-“v=./ 


/— T 


y y 2k 


u!j_ sin 


(A 1 - 26 ; 


^ , TT - / - (Dt sin 4*1 

“lx =/ aTTv " 1 " *>' "v 21 ^ fxj 


/ Ji 
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It be easily aeeb that the Integrand depends on and through 

only the exponential terms which may be integrated by using the following 
w^X Imown integral 


/ e 


-2 

*0},, 


dWy = /n 


(A.1.27) 


Elnall^ we use the following transformation for w and w 

-X lx 


/2k T 


r\ = {b+B cos 


Z = Cb-^’Bj cos 


_ 

— a. 


2k T 


sin p - [b+3j cos 0^) 


] 


f2k T. 


w sin B 

Jl 


1 - f'’*' «/ — jp »lx 


“lx »1 

CAa.28) 


(Di„ sin B 


such that 


\ JCn^+Z^) ; dooj^ « dn dZ 


where the Jocobian of transfoiaation is giyen by 

2k 2k 

J = {(b+B cos ^ sin^ S + (b+R, cos — -s — ^ sin^ B,) (A,! 

<^1 11 1 

As tt® dependence of the integrand on Z is only through an exponential 
tem, the 2 integration can be performed using the standard result 

^2 

'O 

(a, 1.30) 


«> -Z^/J 

/ e dZ = /irJ 


33ms we teve for any f (c^, 1/^ (^Z or ^ 

integral giyen by the 11 fold integral as follows 


29) 


) the collision 
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r K! ,3/2 

^3/2 27TkCTt^+Tt,) 


, 8s ^1/2 


f 8s , 1/2 v2 
3 


/ T|T vjKfif YC“g.n) e3cp {- 


niG 


2kCTt.+'rt.) 


-n^ 


5(2B)^ 
2k Ty. 


5(28^): 

2k Tv7 


} F sin a da dG dn dt. 


(A.1,31 ) 


where 


Ti = r]//j 

g.n = V cos a « n + X 
X =-v(3^ sin cos Pj - B sin({i cos P) 

(Chere are yet 11 fold gf integrations which are to be carried out* 
These can be done only after a particular choice is made for f * 

A.2,1 ^ * C ^ : 


(A.usa) 


Using the Eqns. (A. 1*5) we can write 

r,2 

= - 


(a. 2.1 ) 


m*' 


witere 




s: i.n = — cos a - 

A ^ -m y 


V v’2 


sin <5 sin a 


(A, 2.2) 


“ i-y- " ""lx ' '"x 

Eons. (A.I.II) and U.1.14) are vised in writing down the Eqns. (A.2.2). 

As meGtioned the dependence on 5 is only throu^ 't ^ and the 
integration with respect to 5 is carried out as follows (Use Eqn. (A. 1.15)) 



78 


2tt 


2n {c2} . / {c2} 




(a,2.5/ 


Substituting from Eqns. (a, 2.1 ) and (A. 2 . 2 ) in (A,2 o) we get 


. 2^2 


in. y T 


{- 


3v^-v^ 


ffi 


cos 2a} 


(A. 2.4) 


A. 2.1 Integration with respect to n - 

350 perfOim the integration with respect to n j we substitute 
the value of X from Eqns. (A. 1 . 5 ), (a. 1.9) and (a.I.II) in Eqn. (a, 2,4) 
and in turn substitute this value for } in Eqn, (a. 1.31 ). Ihis 

reduces the n" integration to the following 


/ YCn-A){(gjA^+g2A^)-(2Ag^+3A2g2)’1^Cgl+^Ag2)b^-g2’l^^®'’’^ (A.2.5) 


\^ere 


„ _ o ^ cos 0 ( 

gl=.2 — — .g2 

V ? 


and 


J-, {JL- + cos 2a} 


(A, 2 . 6 ) 


A * V COS ct 4- X 

Ihe above integral gan also be written in the form 


-2 


{ (g^^+g2^)-(2^j^A+3g2^)^ +Cgj+3g2A}Ti^*g2n^}e ^ (A.2.7) 

A 


where 
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g = , 2C ^ ; g2 *= {{Ui) + 

_ _ iT - (a. 2. 8) 

A = A/j ; V = vj/J ; C = (:;~)^ ; rT - 

V 

The integration of (a. 2.7) gives us a function of A, , g^jJ 


3/21 


t2 


” J Cg2 Si A + g2 - j(gl + 3g2 A + gj^ A 


1,- 


t2 


+ g 2 er£ T) 


(A.2.9) 


A^2*2 Integration with respect to a : 

We have from Eqns. (A. 2. 6) and (A, 2.8) 

J a V cos a + X {A,2.10) 

This transformation is used to change the a integration to integration 
with respect to a . In such a ease we have 


sin a da = -dA/v 


V 2 = X -t- V and = X - v 


(A.2.11 ) 


where and are the lower and upper limits (corresponding to 

as 0 and a * TT ) for the new variable A, of integration. 

'Phi « means the following integral has to be evaluated 

3/2 _ 

/ F (A ; gT,g2) dA 

V 

It should be noted at this stage that and g 2 appearing in , 
thaaselves are functions of A, given by Eqns. (A.^.s). fhus afte 
substituting for and ^ 
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3/2 V2 5 _2 

TZ — / ^ I ’ ^-s, \ er£ A} A®- dA 

V Vj, ^.=1 


r3/2 


= -r- OHv 2 ;v,C, 5,^ . H(vj;v,C,C,X)>. 


where 


y^T X? ^ 3 r3£»l -9 

o “ r r ’ ”1 ^ 7T ^ 


p ^ 


p , ag 3X(3g-l) p . 3(5g.l) 7X ^ l-£ , C3g-l)x' 

2 ■ ^ ' ;;;2 ^2 ' 3 ' ^ 2 ^ 2 * ; .2 772 .2 


V- 




2C3?-1)X . ^ 3^-1 

= - — -r— r , rr - — 




?2 ^ 




l-€ ^ (3g^l))^ . 0^ ^ 11 _ X(5U1) 


2v^ 


c — 2 


^2 ' 


5g ~ 1 i , M * C 3g-l)X 

2^2 72 C 


7r2 


2*;^ 2^^ ^ 


^3 


5g.l . 

2-2 ' 


0, = ^i^ 

2^2 v2 


Qr = 0 


n 

•o = 5 


3C35-l)3r2 3(U5) ^ 

_ . + ' '■ .■ < ;■ >»■ •• ~ 

2c' 


^ ; p-i = — ^ " 5 


2 c 2 v 2 


Ro 


. 5(5C^1)X 


R , ^C5g-1} . 2£ ^ H ^ 31^ ^2 
^ 2 r} ^ 


R 2(35.1)X, R, 

^ 2-2 3 2 rfi 


(A.2.12) 


CA,2a3) 
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_ _ 5 

£+1 


n S,+ l 

Po 

. K p 


- h®. Cvj} 


(A.2.14) 


353 e integrals 

X 2 

\.,(x) - / X e dx and (h)^(x) = /* x” erf x dx (A. 2. 15) 

^ 0 

appearing in (A. 2. 14) have been listed in Appendix B 

Integrations in the G space j 

Sqns« (A.1«11) and (A,1,20) give the following transformation for G 


ifeere 


V = G - u, . i 

- - 13 


(A.2J6) 


Bepressating 
vDligae in G 

m 


u. .= u. - u. 
i3 1 3 


G in spherical 


w 

space as 


{A.2.17) 

co-ordinates we have for the elemental 


dG sin dVj dG 

®[X>® Bqn, (A,2,16) y/e can get the following resulfo easily 


(A.2.18) 


v^ » G ^ - 2G u, . cos 7^ + U.2 

IJ 

T » 5 cos Y - U. . 

T - n 


(a.2.19) 
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She flret of these two, gives a transformattoa tor integrating with 
respect to Yj 


sin YidYj- 


G^. 


(A.2.20) 


where G » G/ / J and u • = u ./ /j 

**• J 

Further the Bqns. {A*2,q) and (A. 2. 19) give e 


•p^ 

function of G , 

and 

uT. 

ij 



G 

:: 2 „ 2 

V •" u. . 

2 


u — ir . 

( — _ — ii- ) 

2 V u, . 


(A»2,21 ) 

as the following 


(A.2.22) 


Shus the voltme element d G given by Eqn. (A.2.18) reduce to 


dG * G — — — dG d'f'i (A.2.25) 

u. , 

11 

and the region ot integration is 
: 0 to 2v 

V i ( G - u7 - ) to { G + ) 

11 

G J 0 to » 

Ihe integration with respect to gives a factor of 

region of integration (OBDE) in the 5 , space is ehoMl in Kg- A.Z. 1. 

33iue we have to evaluate the following integral 



(H^ - H^G 




dv 


(A42.24) 



viiiere 


a = n J/2k(T^ + ) 

i 3 


H 2 “ i:2 Cv,C-;c^X) ^ H(v 2 ; v^£,c,x) 
Hi Hi 6r,G;c,^ = HCviiv,C,c,3() 


(a. 2. 25) 


It can be seen from Sqn. (a.2.11), {A. 2.I3), (A,2.22) and (a. 2.25) that 


H 2 £-v,-G;C,X) = H 2 C-v,G; 5 ,X) = KiCv,-G;c,X) = Hi(v,G;C,X) 


(A. 2. 26) 


fhe integral (A, 2. 24) is eqxii valent to the difference of the contributions 
due to the region AHE (Sector of circle with infinite radius and centre 
By bounded by the lines G « 0 and v = G— u^^ ) and the region 
ABC (sector of circle with infinite radius and centre By bounded by- 
lines G = 0 and iT = G + ). Elis is possible only when the integrand 

is well beha-v-ed in the domain ABE (which contains ABC as a sub domain 
in it)y Such that the integral with « limit is not an improper 
integral. careful study of -ihe integrand in (A 2.24 )y that iSy 

the fiinctions gi-ven by Eqn. (A.2.l3)f it can be seen that the integrand 
is bounded due to the exponential damping terms and hence it is well 
behaved, in the domain ABE. 

35ius we can write the integral (A. 2. 24) as 


„ 5^u~ - -7r2 G+u. . 

^ f (H2-Hj^)Ge* dv* - £ dG / CH2-H|)Ge"' dv 

-“ij ° “ij 
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M-D&lly f by effecting the change of order of integrations and using the 
sysmetric properties (A. 2. 26) we get the integral 


r 


v-u,. 





dG 


(A. 2. 27) 


where 


H2=HCv2JV,e,f:>X) 






functions of v , £ , 5 , X and 5 is a function of G, v , 

given respectively by IJins. (a, 2,13) and (A. 2. 22). 1116 integration with 

- 

respect G can be carried out, as is an integrable function 

through the functions which are in turn dependent on | * 

However the analytical integration, with respect to v has been not 
possible. Ihus we get an integral with respect to v as 


/ K (v; a,5,X) dv 


(A. 2. 28) 


where 

H2Cv,a,5,X) - I ^ ^ 


H-l Ji+l 




where 


(A.2.29) 




= / Pg G e dG 


— 2 

%(v,a,?,X) = f 2 G e""^ dG 

_- 2 _ 

I^Cv,a,C,X) == Rj G e’*®'® dG 

V2 - T + V 


(A.2.30) 


’i Y and V + ^ 
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Bie integrals (A.2.30) are listed below 


" 2 


2i;^ 


3Ci X 3X _ C3+2 X)c3 2 ^ 


^2 ‘ c 


^ 2^V 


-2?3X 

P4 =T^ 


; Pc = — 

* 5 2 ~2 

V 


_ ^2 _ El * * E3 




^3(1+^) q ?2 

0 _ 

2^2 ^ 


? 

C 2?2 * ^ 72 


(A.2.31) 


^4 ** _ 2—2 
itrj^ 

?lX 

^ = -i 

O c 


Ro = 


; Qs = 0 

3^3 x2 2^2 
2?2 ^ ^ 
^2-, 


_ 2^1 X _ C3+2x" 3^3 2e^ e2 


V* 


' 


-2^3 X 


R. = 


■4 - ^2 -2 ' S 52 ^ 


i^ere 


g e / ^ e G dG = " ^11 

i n 


52 » C1-5)G e'“®^ df - CEo2 ' ’ ^^12 ' 


■-I 
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6^ 


5j = r {Se-DG -5„)- 


*5(=;) - 2(^ * 5?j) ^Aj (Cp ♦ i Aj(Gn (A.2.32' 

ij a * ‘ 

^op = r »;) and = •^ G, 


a 


!ams v?e have^ iiie collision integral for $ «; c as a seven fold 

3C 

integral given by 


-2nin. 

= iJL { 

ij A 


m 


2k(t^..T,.) 


,3/2 , 8s , 1/2, 8s , 1/2 r ,3 

^kT kT ^ — J 

Kiy. *^1V4 u- . 


ij 


Cv;a,?,X} exp {- |ipi - }K dv d?^ 

Vi V. 


(a.2.33) 


fhle expression I. , (c ) takes o/o form when i = j. However it is 
move oonvinient to carry out the integration in G space* separately 
for this case ( i = j) » 

iioia B^n, (^,2,16 )^ (a. 2.19) and {A.2,22) we have for i « j 


u. = 0 ; V = G ; V = G cos Yi 

11 ^ X ^ 


'"v 2 2 

? = (“) = ^1’ 
V 


(A.2,35) 


Sius in Bi^n, (A.2*12) the functions which are dependent on K 

c^ be integrated with respect to Y ^ to obtain 


CvX?) ' r PfcSin Yi dYj 

Cv*X, 5 ) = f Q£ sin Yi 

(v,X,C) R. sin Yj 

0 


dY, 


(A.2.34) 
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so tbat 1318 collision integral can be written as 


o 2 

T. . fc M 
11 4 
IT 




V. 

1 


s(2E)2 + S(2B )2 

H 2 Cv,a,CjX) V exp {-a } k dv dt. 


2k T 


V. 

1 


(a.2.35) 


wliere 


5 P 

H2 Cv,a,C,X) I -i-i Qf (V 2 ) (A.2.36) 

fc=l i!,+l X ■ < 

are the functions Ejj^ given by (a.2.31 ) 

in which Cj > ?2 ^3 replaced by and respectively 


where 


^ sin dYj = 2/3 


0 


52 = (1-?) sin Yj^ <iYi = 4/3 


0 


(a.2.37) 


53 “ r (35-1) sin Y^ dYj = 0 . 


0 


Kius the collision integral for f - is reduced to a seven fold 

integral both for the cases i ^ i « 3* 

a .3 i= 1/2 j, 


, J_l2 {(b+B cos 4 .)^ (l-Ce.j)^) Ib-B cos +il 


2J 


(1 - (en.n)^}} + I i(P*h * 

- (b + B cos 4»)“ • 


exn 


(A 5 t) 
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jtot&er uslJig the transformation for e , n , w and we can rewrite 
this as 

2 

{f} = <|>j^}^sin^0j+(b+B cos 4))^sin^e] - In (A. 3. 2) 

wbgre I,n,^,p 2^are given by Eqns. (A.I.3 ), (A, 1,8), (a. 1.9), (A.I.IO), 
(A,i, 28) and (A.I.32). In this case {g} is independent of 5 , so 

/ {?> d& = l-n Tfl j 1^] ' i 

0 

Integration with respect to 

Substituting this {|} in Bon. (A.I.31) we find that the Ibllovving 
integral in n appears in i ). 

—2 

1 vCn-'A) * ^ 2 ^ CA-^)^> dn (A.3.3) 


where 


Pi = 


m 


2 ?^ Ji 


2 , 2 

{(b+B, cos dj) sin 


+ (b + B cos <|s)^ sin^ 


B> 



and 


(A.5.4) 


A is given by Bqn- {A. 2, 6) 

ae tategral (A,3.3) can be re«itten with tbe proper limbs for n 
and after earrydng out foe S' integration we get the form 


f{p # . (U2 - 31*1) j? n . t3p, - 2P2) ''' * ^"2 ' “l) " 

A 

ta j5/2 [(3v^-2v2)A+2pjA^J j ' Pi A ] 



,-#5 ’ 


eg 


* * j - 2h2)7: + y a 33 erf A) 


* (A ; yj, y^) . 

lateggatloia with respect to « 

■Jae saae traasforasation used for- liie earlier case, given by Eqti, (A. 2.1 
is adopted for the ce integration, to obtain the integral as 


-r~ / = ~-rr' ’ (A.5.6) 


V V. 


where 


Sp =■ S(v ; Uj.Uj) = f ^ [OMj - 2^2) Vp + Uj vd] * 


ij- {4pj . 3ui) - 4(3wj . 21,2 ) - 4mj v^] CTf [14^2 ' 

-vj (A.3.7) 

- — vie 
8 P-* 

Integration in G space : 

Here again the integrations are done by adopting the transformations 
used in the earlier case given by Eqns. (a. 2.19) and (A. 2. 23)* Ihus 
we hs.Te find the integral 


,7/2 


u< . 0 

ij 


d€ 


CS2 - Sj^) G e dv 


(A.3.8) 


u 


13 


ttese integrations can be done analytically, wi-th out effectl g 


the change of order of integrations, unlike in t/he case 
with respeot to v we get from (a.3.s) 


of ? = 


5pi)Vp 
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\i. . 

13 


*c& 

/ CSj - P^ ^ Sj - S^) ff o- 




dG 


(a.3.9) 


5Ba«r@ 


?Kp : Bi-”2> = «ric - 2P2)fi^ •* 3Pj S®T 


1 ^ I 

p- 


1 


* 240 - 2ii2)^p - i2y G^] (A.3,10) 


1 2 
erf ^ ff40i£_ « 24is.'i + rAn,.- _ \7T^ n.. 7^^ _-G 

F 4S0 


[(40^2 - 24p^) + (40^2 - 54y,3G^. I2y G^] e'^P] > 

P * P 


Ej - G 4. X^ } ^2 ~ ^ ^2 

SI « - 5” + X~: S’, = - 5^ + K, . 


(A. 3. 11) 


Ifitegra'^a with r&spect to G is carriei out using the transformations 
»5.1l} azi^ after rearranging the tejrjs we get from (A. 3 . 9 ) the required 
iSl SB 


J/2 


II* - 

13 


{S(Xp;a,C} - ? (X. ; u,^) 


(A.3.12) 




s = sex I = ~ {^£- 
F p * ^2 120 


[5(c-3){x 4 Ca,Xp) 


- Xj XsCa.X^)} + 3(5-1) 1X6 («.y ' ^ xC«,Xp)}J 
"■ 240 (5'2f“'V - ’i'I<“>V* 


♦ 20(5- 3) (Xj, (--Ca.Xp) - '»4(a,Xp)l + 


ff‘ 


{X„ (a,X^) *6(a.Xp))] * 240 ^ 



- (27-7?) 

it^ + 6CC-l){X^2^Ca,X^)-Z^(a-Xp)}]> 
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ae int©gpale 


5^C«»x) 


^-OCG-x}" ^n 


6' 


m 


C«,x> 


/“ g” e-«' dG 


{A,3.U) 


as:^ 


XCafX} 


/“ G" 


erf G 


dG 


baT0 eval'oateS and glrsJi in Appendix B, 

b^afore arriving at (a, 3.12), the quantities and P2 appearing 
in are e3q>2^en3©d in tems 01 C using Dqns, (A.I.IO) and (A. 3*4)* 

Ssleg -gie issmilt (A.?, 12) in Ebn. (-'.I.3I) we get the coliisien integral 
for HilB case ss 


h} % *^1 = 


X 


X f. 




■2k(T 

J 


{ 85 j.1/2 |- Ss ^ 1/2 b' 


kT„ 


kT 


u 


13 




CK' { - 


s(2B)^ 

2kTy, 


2k T^_ 
j 


2 

-} df^ 


(A. 3.15 


a six fold integral which has to he evaluated numerically. 
m In c ^ casej. ?&en i = 3^ this reduces to a o/O form. 

'Si s colli sloa integrals for such a case v/e follow^ 


2 

ae followed lor the c case. 

.A 
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^ of , the integration with 

ki "Vj gi^e a f&o^or 


r 


siit V Ti ~ 


«e gat the i^ilisic-i: integral 




11 ' 1/2 > 1-2 r“ ,2 


r. X o-/ t, . o^> V , z f y 

^kT’ ^ " «{''2>Vi,P2) 


‘t • 
1 


'v 63 ^ {-iv^ . 


s(2B)^ -- silZ. y 

— — } K dv dt, 

2 k Ty. 6 


t^feich IjB oh^ously a seven fold integral. It has not been possible to 
o&XT^ut the V integration analytically for i = j case. 

A .4 I * 1/2 3 ( 2 B)^ : 


,(A*t, 5 ) '»e oan obtain for this case 


{f } = - I V (Ij sin cos 6 j - B sin ^ cos B) (A. 4.1 ) 

As this is Indep^adent of *5 in this case we get 

IfF = {.t} 

^t<^raMoa wlHi respect to n ; 

for I frcn Ifen, (A.1.8),(j^.1.9)j (A,1.10)j(A. 1.28) and 
(a* 1*32) get ■l^e follo’ffing integral in n 


H 


V (A-y)*' e 


A 



dp 



¥ 



(I? sin t cos B 


sin cos B) 





X 


(A.4.2) 


(A.4.3) 
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lategratijag (A.4.2) we get 

^ ( 1 + 2 !^) - ~ (1+2A^) erf A - ^ Ae" 

Integration v?ith respect to a ; 

Sis is done by using new variable A, as indicated earlieor 
A integration we get 


Jv ► 


V 


where 


Vj ^ 


V 

P 




- \ 0 


-/ 

- v^) - 


^ (3Vp 2v^^3 erf v^} 

Integration in G space ! 

Here again the procedure used in the earlier case is adopted 
the required integral as 


,5/2 — 00 u+Uj ^ „ 

~r~^/ ^^2 - V ^ ® 

«ij 0 

Integration with respect v gives us the above integral as 


•“2 


ail. /" [V, - V, * Vj - V^)G 6-“®' dff 


u. . n 

13 


} {A. 4 . 4 ) 


After the 


(A.4*5) 


(A.4»6 ) 


to obtain 


( 1 * 4 . 7 ) 


(A.4.8) 


where 
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*i) erf G^, 

and are same as given by Eqn, («-,3.1l), 

Cariying out the integration and rearranging teims gives the 
Integral (ii.4.8) as 


C^. 4 . 9 ) 


3/2 - 

J V 


u. . 
13 


{VCa»X2) - V(a,X,)} 


where 


(4.4.10) 


% = V(a,Xp) = [(Xjff.Xp) -5Xp X2(«,Xp)} 


{Xj Cs.X^) - X^Cti.X 1}] 




-Xp -PjCc.Xp)}] - [5{Zj(a,X^) - X Z„(a.X n 


+ 12{23Ca,Xp3-Xp Z2t«»V^ ^ ^ ^ 


S’ih^lly we can write dovm the collision integral for this case as 


,2 -2n4 n- 


I ^ 3 f 

^9 2 ^4 2 k(Tt 3 Tt.) 


,3/2 r 8s ,1/2 , 8s 1/2 b 




/ 7 (Vj-Vj) exp {- |:2|L - - -■/ ■■■-> ih 


2k Tv 2k Ty. 
1 J 


(A.4.12) 



APPEtJDIK B 


Some of the important Integrals used in Appendix A have been listed 
here. 


_ 2 

3.1 Integrals = /x^ e ^ dx ; 

Aq = erf X 

Aj " ^ 2 

A = (erf X - X e"^ )/2 
2 

1 / 2 \ ~x^ 

= - t (1 + X ) e 

2 , 

= t{ 3 erf X - e"’' (3x + 2x^) > 

2 

= -■^(x4-2x + 2)e 

2 

A^ = ■! {15 erf X - (I5x + 10x^+ 4x^)e“^ } 

6 o 

B.2 Integrals #^(x) = / x^ erf x dx : 


= X erf X 4- e ^ /2 


=! i {(2x'^-1 ) erf X + xe**^ ) 


= ■“ {(2x^ erf X + (x +l) e > 

6 

-x^ 

(4x^-3) erf X + (2x^+3x)e ^ > 

1o 

= n{2x5 erf X + (^2%^+ 2 ) e"" > 


— {(8x^- 15 ) erf X + (4x^+ 1CK + 15x) e > 
48 


5 
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B*3 Integrals 

X„c«,x) = r g" dG 

»= /rtTa 

5(2= (1 + 2- x^) X^/2a 

X<r = (3 + 2 ct x'^) X X /2 Ct 

^ o 

= (3 + 12a x^-»- 4aSc^) Xq/4«^ 

— (15 + 2Gn X -f 4^^x^)x Xq/4 

Xg = (15 + 90 a x^+ 60 + 8 o>\^) X /ea^ 

0 

B.4 Integrals i>^Ca,x) = /” G^ 

^ca 

= »^7Cl+a) exp {-c!X^/(l+a)| 

4>i = ax 4^g/(l+a) 

^2 ~ Cl+«+2a^x^) ^Q/2(a+l)^ 

" {3a (l+a) + 2a^ x^} x ii/^/2(,l+a)^ 

^4 = {3(a+l)2 + 12 (a+l)a^ x (a+l)"^ 

^5 {ISaCi+a)^ + 20a^ (l+a}x^ + 4a^ x"^} x ^Q/4Ca+l)^ 
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B,5 Integrals Z:(a,x} = erf r 


^ puttlce y = )?■ fonortng results are listed. 


X eri y 
o ^ ^ 


= (e’^ /2 v'SX»TT“ I + X erf y) 

= {iii±i!2 erfye. i?? .J- L^. e'^ 


2 t^oi(l + 


e } 


r (2a X + 3)x 
1 1 :LL^ erf y + 

2a 

4 - { 3 ’f 6 x^)a^+ (5 + 2 x^)a+ 2 } -y^ 

4 (a 4- 

r / 2 4 2 » 2 

t v4 c*"x + 12a z +5) erf y/4a 


J Xo 


+ f ax\ o!* 4 (l2x\ 12x) a^+ (ax^+ 26x) a^ 

2 

+ ( 2 x ^4 13 x)a + 5x } /{ 4a^/^(a + 1 >] 

O A 0 2 

I { 4« X + 20a % 4 15} X erf yA^ 

+ C (15«^+ 20a -*■ 8) (1 + + 2x^ A + 0 “ 

_ 2 

( 50 / 4 5Ca^ -f 35a + 9) + 4x^«^ (5«^+ I 0 A lOa + 5^ +1 )}e ^ ]x^ 

(1 erf y + 3e ^ ) Xq /(“ + 1 


^iiere 





A » 




3 -It ^ 2 4-t? 2 2 ^ 

.|. a ^1 + 2} rt -f a + ct^ (a -{- 1 ) 


15 „2 


a'"^ (a + 1 .f ^ (a + l) 


^ 8 ^ 
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^ 45 2 ,.,. 2 / ^ 5 2f ^2 5 , 

» {-r‘‘2X + 1t?a ('i + 1; s + 3 a (a + 1 ) x } 

»^(a +'1' ) - -f -r x^ (a -f- 1 ) 

O C. 

+ X^ a^(l •+ a)^} y -t- {15xot+ lOx^ 

a(a 4- 1 )} (l + y^)/ a (a + 1 ) - { -^ a + 

»-p i 2 

«(s + t}} (3y + 2y^) + { 3x (y^ + 2y + 2) 

'i a(a 'i t ) - 4 (15 4- 10y^+ 4y^) y * 
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appmbix a 


C .1 Collision Integrals D. . and U t 

- 13 0I-I 


Eroia Chapter 5 we have 


13 


= /rf.,fpd't 
= / f.f, r{- g.a) K a a dt, dt 

i J . «> ^ A 


(c,l*1 ) 


yellowing the procedure given in Appendix A, 16 integrations were 
carried out for the case i j,, to get the six fold integral 


hi iij 


in 


4 ^^ 2 x!:(T^.+T ^.3 


J-./ {R^Ca,Xj,)-Rs(a,X2)}e)<p { 


s(2B) 


2 sC2^i)' 


2k T, 


2k 7 


-} sin g sin By dB dpj^ dB dBj^ d^ di>y 




(C.1.2) 


where 

R^Ca.X) = ^ 2" ^ X3) ' ^ h 

* 6 Zj - 4x Zj » 4 Z 4 ) - iy t-2x 

. 2 !(/j - 2 i !<2 + 21)13) (0.1 .5) 

sse q.uactttles u. . S,X,,Xj, and J appearing In Eqn. (O.a.2) nave 
already been defioed in Appendix A and the integrals x„ , ^a.'^n 

listed in Appendix B- 





if'' 


Jbr 1 = 3 we have 
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D, 


= ^nf 2 }V2 r 

>1 « * 


4irk T* 


D v^ ■ s(2:i)2 + s(2B )2 

^sd — - -I — .]} 


2k T, 


sin B sin dg dS^ dB dE^ do 


(C.1.4) 


•vbe3?e 


ift 




(>x,x) . ^ CX3 - X X 2 ) - J Pj - 2x ZjPj-xZ^^i (0.1.5) 


Here a^in all the quantities used are defined in either Appendix A 
or B» Hie collision integral in the atoiaic species conservation 
equation may be obtained as a special case of D. . by effecting following 

J- J 

substitution as indicated in Ci^pter 2 or 3* 


= 0, - 0, to = 0, g = 0 


(G.1,6) 


0*2 Collision Integrals fh -(§_ ) and ^ ^ * 

2 

'.fe can write the expression for integrals s{2B) ) straight 

away using results of section C.l as follows* 

1 / j 


-r 


2 H 


m 


^ 4 vV 2vk(Tt +T^ ) 


3/2 1 


u. 


't, 
1 T 


IJ 


/ CR_C«,XJ - K(a,x^) dsczroh exp {, 


2k T. 




s{2B f 

sin B sin dB dBj dE dBj d(ti d<f)j 


CC,2*1) 
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agi fer 1 = 5 

c|st2nS= - n; j 

4nlcT^jl^ 

^Sd %s( 25)2 j ex^ {[ — --- ■■ i-]> 

- 

sin a sm da dgj^ dB dB^ di^ d<;.| {c, 2 . 2 ) 

2 

Box the integrals R (? Jh “j_) the method indicated in Appendix A 
is foUwed aicd ,ve can write for the final six fold integral for the 
case 1 ^ 3 as 


Rij (ijiuf) - I (ij, ..2) dt 

^3/2 _i_ j j3/2 


"l 


“i/ir 2irk(T^ +T^, ) "13 

i j 

{Py {ci,Xj) - (5,X2)} exp {- 


s( 2B)2 

2fcT„ 


where 


s( 2 B,)^ 

} Sin g sin dg dg^ dB dE, dij dijj^ 


P-j.Ca,)f) = 3 (X4 - X X3) - [Zj - C3-x) Z2 


(C.2.3) 


CC.2.43 


(X + i) Zi + f ZJ 1 [ I 'I'o + ^1 - *o]> 


Bir i = 3 we have 


¥*■' 

m. 
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J, ...^) = . nj ( !D_)3/2 J y 


s(2r)^ + 5(2tj)2 

r,x) (*• — . — 1 sin ^ sin e, 

2k T„ ' ^ 


(C.2.5) 


d8| dU <53 J d^v dfj< 


i^ict tk^ integrals 3 . .^(5) fiXJia S. .(< 5 ) the same substitution 

Oi2 13 - 

in aqn« (C*1 .6 ) is used. 

have all the required collision integrals for the 


iaEliag case. 
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